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1.  INTRODUCTION 


The  work  reported  here  is  concerned  with  the  use  of  the  cepstrum 
for  the  reduction  of  multipath  interference  in  communication  systems.  It 
is  basically  a continuation  of  work  performed  earlier  and  reported  in  [1]. 
The  application  toward  which  this  effort  is  directed  is  HF  data  communi- 
cation. Through  the  use  of  the  cepstrum  it  should  be  possible,  under 
suitable  conditions,  to  extract  at  the  receiver  a good  estimate  of  the 
pertinent  multipath  information  on  a continuing  basis,  without  recourse 
to  interspersing  a special  channel  sounding  signal  with  the  normally  trans- 
mitted communication  signal.  This  multipath  information  could  then  be 
used  to  adjust  a filter  which  is  designed  to  cancel  the  multipath  inter- 
ference in  the  received  signal. 

The  earlier  work  involved  the  development  of  suitable  programs  on 
a large  computer  - the  GE  645  (now  Honeywell  645)  - and  thtir  evaluation 
using  synthesized  signal  sample  sequences.  The  results  which  were  obtainea 
looked  promising.  The  present  effort  was  therefore  undertaken,  with  the 
aim  of  developing  and  evaluating  on-line  cepstrum  computation  using  a small 
computer  - the  PDP-8.  Such  an  arrangement  allows  the  processing  of  actual 
modem  signals  that  have  been  subjected  to  live  or  artificial^ imposed 
multipath  distortion.  It  also  differs  from  the  earlier  work  by  bringing 
into  play  the  effects  of  very  limited  word  length,  the  effects  of  sampling 
distortion,  and  lack  of  synchronization  between  the  signal  source,  the 
multipath  delay,  and  the  sampling  process. 

The  performance  of  such  a cepstrum  processing  scheme  depends  to  a 
certain  extent  on  the  waveshapes  that  are  transmitted.  Whereas  the  pro- 
cessor can  be  thought  of  as  strictly  supplementary  to  an  existing  communi- 
cation system,  the  possibility  of  a complete  system  design,  requiring  the 
specification  of  optimum  waveforms,  was  also  considered. 

For  greatestease  of  presentation,  the  main  body  of  the  report  is 
rather  brief  and  in  summary  form,  while  the  technical  details  are  presented 
in  a number  of  appendices.  The  basic  concepts  of  cepstrum  analysis,  and 
its  use  in  processing  multipath  signals,  are  sketched  in  Section  2,  and 
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Section  3 contains  a brief  review  of  the  earlier  work.  The  current  effort 
and  the  results  obtained  are  discussed  in  Section  4,  with  reference  to  the 
Appendices.  The  principal  conclusions  are  listed  in  Section  5. 
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2.  CEPSTRUM  ANALYSIS  APPLIED  TO  MULTIPATH  SIJKALS 
2.1  The  Complex  Cepstrum  - Continuous  Case. 

A cepstrum  of  a signal  x(C)  is  obtained  upon  applying  the  following  se- 
quence of  generic  operations: 

Fourier  transform  - logaritta  - inverse  Fourier  transform  (2-1) 
Interest  In  the  cepstrum  Is  due  to  the  fact  that  If  two  signals  /:^(t),  *2^) 
are  convolved  In  the  time  domain,  then  under  suitable  conditions  tbs  cepstrum 
of  this  convolution  consists  of  the  simple  addition  of  the  cepstz-*  u the 
Individual  signals. 

The  detailed  specification  of  (2-1)  can  take  various  forms,  leading  to 
different  versions  of  cepstra.  One  such  version  has  been  called  the  "complex 
cepstrum"  [4]  because  the  complex  logarithm  is  used  in  (2-1),  although  the 
cepstrum  Is  actually  a real  function  (If  the  signal  is  real) . Given  a signal 
x(t),  then 

QV 

Its  Fourier  transform  is  ^(f)  * j x(t)e”J^ir^tdt;  (2-2a) 

Its  logarithmic  spectrum  is  £ (f)  - cln9r(f)  (cln  '*  complex  logarithm); 

x x (2-2b) 

m 

its  complex  cepstrum  is  €^(t')  » | •56j(f)e^^xft  df.  (2-2c) 

The  multi-valuedness  of  the  complex  logarithm  is  resolved  by  requiring  con- 
tinuity of  £ (f).  The  signal  x(t)  is  recovered  from  6?  (t')  by  applying  the 
x _i  x 

sequt  ce  of  operations  F.T.,  exp.,  F.T. 

Now  suppose  the  signal  x(t)  is  transmitted  over  a linear,  time-invariant 
channel  with  Impulse  response  h(t) . The  received  signal  is  then  y(t)  ■ 
x(t)*h(t)  and  we  have 

tFy(f)  - ^(f)  ^(f),  (2-3a) 

^.(f)  - (f)  + ^(f)  (2-3b) 

and  also 

Cy(t')  - Cx(t')  + ^(t’)  (2-3c) 

Thus,  In  both  the  logarithmic  spectrum  and  the  cepstrua,  the  contributions  of 
the  transmitted  signal  and  of  the  channel  add.  If  the  impulse  response  h(t) 
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is  Impulsive,  then  so  Is  C^(t'),  so  that  the  characteristics  of  the  channel 
are  most  easily  identified  in  the  cepstrum  domain  (assuming  that  the  signal 
cepstrum  is  not  also  impulsive). 

Two  different  ways  in  which  the  cepstrum  can  be  used  to  reduce  multipath 
distortion  are  (a)  "homomorphic  filtering",  and  (b)  channel  identification 
coupled  with  a cancellation  filter.  In  the  firnt  case,  linear  operations  are 
applied  in  the  cepstrum  domain,  in  order  to  reduce  - and  possibly,  eliminate 
altogether  - the  effect  of  <2^(t')  in  (2-3c).  Transforming  back  to  the  time 
domain  results  in  recovery  of  x(t)  with  little  or  no  distortion  due  to  the 
channel.  In  the  other  case,  the  pertinent  channel  parameters  are  estimated 
from  6?y(t')  and  these  estimates  can  then  be  used  to  adjust  a multipath  can- 
cellation filter. 

The  simplest  type  of  multipath  channel  is  one  which  introduces  a single 
discrete  multipath  component,  with  amplitude  and  delay  that  is  approximately 
constant  over  an  observation  interval.  The  received  signal  y(t)  is  then 

y (t)  - x(t)*m(t)  - x(t)  + ax(t  - t),  (2-4) 

where  the  multipath  channel  impulse  response  is  m(t)  » 6(t)  + a5(t  - r) . 
Applying  (2-3)  to  the  received  signal  (2-4)  gives 

%(f)  « Tx(f ) (1  + ae~j2nfr)i 

y A 

-J27rfT, 


£y(f)  - (f)  + cln(l  + oe‘J"lT) 


- ^(f)  + ae-j2*fT  - 4 e"j47rfT  + f-  - ...  (for  a < 1) ; 

^ 2 3 

. e tf(f)  - e (t*)  + o6(t'  - t)  - f-  6(t'  - 2t)  + 6(t*  - 3t)  - ... 

y 1 J (2-5) 

The  effect  of  the  multipath  in  the  complex  cepstrum  is  a sequence  of  impulses 
with  spacing  t and  weights  determined  by  a.  Therefore,  if  6*  (t‘)  contains 
no  impulses,  the  multipath  parameters  a and  t are  easily  identifiable  in  the 
cepstrum. 

2.2  The  Complex  Cepstrum  - Discrete  Case 

In  practice,  digital  processing  is  employed.  Given  a sequence  of  N 
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samples  {x^},  k • 0,  1,  . ...  N-l,  then  Its  discrete  Fourier  transform  (DFT) 
is 


(dx(n)}  - 


N-l 

!i  r 
In  l 
k-0 


xke 


-J2^' 


n ■ 0,  1, 


, N-l; 


Its  (discrete)  log  spectrum  is  (tx(n))  - cln{dx(n)); 


its  (discrete)  complex  cepstrum  is 
/N-l 


(«Mk)} 


l * (n)  e 


k * 0,  1, 


N-l. 


k-0 


The  multl-vsluedness  of  cln  is,  in  general,  not  readily  resolved  in  a unique 
manner  in  this  case;  however,  this  is  of  no  concur^  at  this  point.  Suppose 
y(t)  in  (2-4)  is  sampled  N times  over  0 < t < T,  giving  the  sample  sequence 
(yk),  k-0,  1,  N-l.  Then  (y^)  - (x^)  + a{xk_r>,  where  negative  sub- 

scripts refer  to  samples  taken  prior  to  t - 0. 


At  this  point,  in  order  to  simplify  the  discussion,  the  following  ideal- 
ized case  will  be  assumed.  Namely,  assume  that  t - rT/N,  where  r is  an 
Integer,  |r|  <<  K,  Also,  assume  x(t)  is  periodic  with  period  T,  or  at  least 
x(kT/N)  - x((N  + k)T/N)  for  -r  < k <.  -1. 


In  this  idealized  case,  transformation  of  tbs  sequence  {y^}  gives 


-j2' 


nr 

r"N, 


{d  (n)}  - (d  (n)(l  + oe  ")} 

y * 


(2-6a) 


U (n)}  - U (n)}  + {cln(l  ae 

y *• 

nr 


-32' 


nr 

Tr. 


)) 


"J2TriiN  a2  “j4lrDN 

- Ux(n)}  + {ae  N - f-  e N ♦ ...  } 


(2-6b) 


(Cy(k)}  - {cx<W)  + {cm(k)J; 

where  {cn(k)}  is  the  cepstrum  of  the  multipath  channel  response. 


(2-6c) 


c (k)  - 

m 


- Nl--g)!. 

p • 

0,  otherwise 


for  k - pr  mod  N (p  - 1,  2,  ...) 


(2-7) 
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Thus,  for  the  idealized  case,  the  discrete  cepstrum  (2-6),  (2-7),  exhibits 
the  same  features  as  (2-5)  for  the  continuous  case  (Fig.  1).  The  impulses 


Fig.  1.  Portion  of  Discrete  Multipath  Cepstrum,  for  N ■ 256, 
r » 8 , a » 0.5 

in  (2-5)  have  become  "discrete  impulses"  in  (2-7).  (The  effect  of  removing 
the  idealizing  assumptions  is  considered  in  Section  2.4.) 

Extractabillty  of  multipath  Information  from  (2-6)  no  longer  depends  - as 
in  (2-5)  - on  the  absence  of  impulses  in  the  cepstrum  of  the  desired  signal; 
rather,  the  variability  of  (cx(k)>  must  be  small  over  that  range  of  k in  which 
the  most  significant  component(s)  of  (2-7)  are  expected  to  fall.  If  this  is 
the  case,  a good  estimate  of  a and  the  value  of  r can  be  obtained.  This 
allows  an  approximation  of  {cn(k>}  to  be  subtracted  from  (cy(k)},  leaving 
(cx(k)}  with  little  distortion.  Inversion  of  the  cepstrum  computation  then 
yields  {x^}  with  little  distortion.  Alternately  the  estimates  of  a and  r can 
be  used  to  adjust  a multipath  cancellation  filter,  in  which  case  it  would  nrt 
be  necessary  to  perform  a cepstrum  computation  for  every  consecutive  sequence 
of  N signal  samples. 

2.3  The  Real  Cepstrum 

It  turns  out  that  for  typical  communication  baseband  signals,  the  vari- 
ability of  (c  (t)}  is  not  small  compared  to  {c  (k)}.  so  that  a useful  esti- 

X — — CJ 

mate  of  the  multipath  parameters  a and  t is  not  obtainable  from  the  complex 
cepstrum.  However,  excessive  variability  of  {cx (k) } is  only  exhibited  by 
the  odd  part  of  the  sequence,  which  is  obtained  from  the  imaginary  part  of 
(1  (n)}  (the  phase  spectrum).  If  the  odd  part  of  {c  (k)}  is  deleted,  there 

* y 
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remains  &r{Cy(k)},  the  even  part,  which  Is  obtained  from  (&t{i^(n)}  and  Is 

therefore  called  the  "real  cepstrum".  From  (2-7)  It  can  be  seen  that  the 

real  cepstrum  of  the  channel  response,  Bv'lc  (k)},  incorporates  all  the  multi- 

m 

path  information  except  the  polarity  of  r. 

For  these  reasons,  only  the  computation  of  the  real  cepstrum  has  been 
Implemented  as  part  of  the  work  described  in  this  repor*.  The  notation  of 
(2-6c)  is,  therefore  used  from  here  on  to  denote  real  cepstra,  thus, 

Ux(n)}  » ln|{dx(n)}|,  etc. 

The  ambiguity  in  the  sign  of  r,  resp.  t,  which  is  inherent  in  a multi- 
path  estimate  based  solely  on  the  real  cepstrum  is  assumed  resolvable  by 
other  means.  The  basic  processing  steps  are  therefore  as  shown  in  Fig.  2. 


Fig.  2.  Computation  of  the  Real  Cepstrum. 

Since  the  real  cepstrum  does  not  carry  the  complete  signal  information, 
it  cannot  be  used  to  reconstruct  the  undistorted  signal.  Therefore,  the  pri- 
mary effort  of  this  investigation  has  been  to  evaluate  the  extractablllty  of 
the  multipath  parameters  (o  and  t)  from  the  real  cepstrum  of  the  received 
signal  1 

(c  (k)}  - {c  (k)>  + (c  <k)>.  (2-8) 

y x m i 

The  following  observations  can  be  made  with  regard  to  the  extraction  of  a 

and  t : i 

a)  Since  {c  *(k))  is  even,  only  the  first  N/2  cepstrum  points  need  to  be 

examined . 1 

b)  If  the  range  of  possible  multipath  delays  t is  known,  the  relevant  por- 
tion of  (Cy(k)}  may  be  further  reduced,  since  only  the  first  spike  (or 
the  first  two  spikes)  of  (cn(k)}  are  generally  sufficiently  large  to 
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contribute  to  detection. 


• c)  {c  (k)}  is  a random  sequence,  since  it  depends  upon  the  transmitted  in- 


formation. However,  a mean  cepstrum,  {cx(k)},  for  a given  type  of  trans- 


mission, can  be  obtained  from  measurements  and  can  be  subtracted  from 


(Cy(k)}.  This  leaves  the  various  cepstrum  points  as  approximately  zero- 
mean  random  variables  (except  for  the  superimposed  multipath  cepstrum). 


d)  The  first  spike  of  the  multipath  pattern  (p  * 1 in  (2-7)) is  then  identi- 
fied by  locating  the  largest  magnitude  cepstrum  value  within  the  relevant 
region  (after  subtraction  of  the  mean  cepstrum).  Either  a single  spike 
can  be  used  to  estimate  a (in  which  case,  the  value  of  the  spike  divided 
by  N/2  is  the  estimate),  or  a composite  estimate  based  on  the  first  two 
spikes  can  be  obtained.  Only  the  first  scheme  has  actu.-lly  been  imple- 


mented. 


e)  A detection  threshold  is  needed  which,  although  causing  some  cases  of 

weak  multipath  to  be  missed,  minimizes  the  chance  of  a false-alarm.  Espe- 
cially, the  chance  of  combined  occurrence  of  a miss  and  a false  alarm 
(at  an  incorrect  value  of  t)  should  be  minimized. 


2.4  Arbitrary  Signals 


Departing  now  from  the  idealized  assumptions  introduced  in  Section  2.2, 
suppose  first  that  x + r(T/N);  i.e.,  x is  not  an  integer  number  of  sampling 


intervals.  It  can  be  shown  [1)  that  in  this  case,  {cm(k)}  is  approximately 


expressed  by  samples  of  sinc-f unctions  centered  at  x,  2x,  ...  (Note:  sine  t = 
(sin  irt)/rt) . Thus, 


Cm(k)  « 


OO 

-N  l sincfeT  f - k)mod  N)  • 

■ p«l 


(2-9) 


In  order  to  extract  the  multipath  parameters,  it  is  therefore  necessary  to 
interpolate  on  sinc-functions.  The  basic  steps  comprising  the  detection  al- 
gorithm for  extracting  estimates  oi  and  x from  the  largest  multipath  cepstrum 
spike  are  then  as  follows: 


i)  Over  the  relevant  portion  of  the  cepstrum  (see  item  b above)  subtract 
the  mean  cepstrum. 

ii)  Identify  the  largest  cepstrum  component  that  exceeds  threshold.  If 
none  is  found,  a "no  multipath"  decision  is  reached.  Otherwise: 
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Hi)  Determine  whether  one  of  the  cepstrum  points  adjacent  to  the  largest 

A 

multipath  component  also  exceeds  threshold.  If  not,  and  c(k)  denotes 
the  largest  magnitude  cepstrum  component,  then  a ■ 2c(k)/N;  t - £t/N. 
Otherwise: 

A A A 

iv)  Interpolate  Detween  c(k)  and  either  c(k  + 1)  or  c(k  - 1),  whichever  is 

AAA  A 

larger.  Let  $ ■ c (k  + l)/[cy(k)  + cy(k  +1)1;  then  a » 2cy(k)/(N  sinc$) , 
and  t - (k  - 1 + <fr)T/N. 

Practical  communication  signals,  of  course,  do  not  have  the  periodic 
property  postulated  as  part  of  the  "idealized  case"  in  Sec.  2.2.  If  a re- 
ceived signal  (2-4)  is  observed  over  a finite  interval  (0,T),  it  is  therefore 
no  longer  of  the  form  x(t)  + ax(t  - t).  Only  the  portion  of  the  signal  x(t) 
within  (0,T-t)  appears  delayed  (Fig.  3).  On  the  other  hand,  a portion  of 
x(t)  outside  the  observation  interval  enters  as  a result  of  the  delay T into 
the  slot  (0, t) . 


Fig.  3.  A Signal  and  Its  Delayed  Version 
Within  an  Observation  Interval. 

The  effect  of  this  phenomenon  on  the  cepstrum  can  be  inferred  in  the 
following  way.  Writing 


fx(t), 

0 < t < T - T 

xl(t)  - 1 

lo. 

otherwise 

(x(t)  , 

T - t < t < T 

x2(t)  " { 

lo. 

otherwise 

f x(t) 

1 

H 

A 

ft 

A 

O 

3 lo. 

otherwise. 

( 


the  following  expression  is  obtained  for  the  received  signal  in  (0,T) : 

y(t)  » x^t)  + ax^(t  - t)  + x2(t)  + ax3(t  - t).  (2-10) 

Only  x^(t)  carries  the  multipath  information  (in  the  form  of  (2-4)).  The  re- 
maining terms  appearing  in  (2-10),  x2(t)  + ax3^  ~ T)»  can  therefore  be  thought 
of  as  additive  "noise"  superimposed  on  the  transmitted  signal  x^(t).  Any  ad- 
ditive noise  in  the  time  domain  becomes  combined  non-linearly  with  the  signal 
cepstrum  to  produce  Increased  variability  in  the  cepstrum,  and  it  also  de- 
presses the  magnitudes  of  the  multipath  spikes. 

The  latter  effect  can  be  understood  by  considering  the  effect  of  an  addi- 
tive disturbance  c(t)  in  the  DFT  domain.  At  a frequency  point  n where 
|dx(n)(l  + ae  J^irnr/Njj  £or  “idealized  case")  is  very  small, 

addition  of  a small  random  vector  d£(n)  is  more  likely  to  increase,  rather 
than  decrease,  the  magnitude  of  the  DFT.  This  Increase  becomes  magnified  by 
the  logarithm.  Since  the  valleys  of  {jl  + ae~^irnr^N|)  in  general  correspond 
to  smaller  values  of  {|dx(n)(l  + ae  J^Trnr/Njjj  tjian  the  peafcg,  it  can  be  seen 
that  the  additive  component  (d^(n)}  tends  to  lower  the  estimate  of  the  multi- 
path  amplitude,  a. 

Since  this  disturbance  is  introduced  particularly  at  those  frequency 
points  where  d^(n)  is  small,  a choice  of  sampling  rate  is  motivated  which  is 
not  too  high,  so  that  the  DFT  domain  does  not  cover  an  extensive  frequency 
range  containing  little  energy. 

For  a given  sampling  rate,  the  extent  of  the  disturbance  produced  by 
these  effects  in  the  cepstrum  is  minimized  if  t is  as  small  as  possible  com- 
pared to  N.  This  is  one  motivation  for  using  a large  N in  computing  the 

cepstrum.  It  has  also  been  possible  to  considerably  reduce  this  Interference 

2 

through  the  use  of  a data  window  with  a sin  -rise  ("Hanning")  over  the  first 
and  last  few  points  of  the  sample  sequence. 
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3.  SUMMARY  OF  EARLIER  WORK  [1] 


A computer  study  to  assess  feasibility  of  the  processing  scheme  de- 
scribed above  was  carried  out  on  the  Honeywell  635  computer  at  RADC.  For 
this  purpose,  various  versions  of  cepstrum  computation  and  multipath  extrac- 
tion programs  were  written  and  applied  co  synthetic  signal  sample  sequences. 
The  signal  sample  sequences,  which  were  generated  by  computer,  simulated  the 
sampled  output  of  a 16-channel  HF  modem  in  accordance  with  MIL-STD-188C . The 
output  of  such  a modem  consists  of  16  frequency-division  channels  of  75  baud- 
per-second  quadraphase  PSK,  with  center  frequencies  equally  spaced  from  935 
Hz  to  2585  Hz.  The  computer-generated  signal  sequences  simulated  independent 
random  modulation  on  all  channels,  and  were  based  on  various  sampling  rates, 
however,  8175  saraples/sec  was  used  for  most  tests.  Sequences  of  128,  256,  512, 
1024,  2048,  and  4096  samples  were  used  (resulting  in  cepstra  of  the  same 
lengths),  but  most  of  the  work  was  done  with  512  samples. 

A single,  discrete  multipath  component  was  introduced  into  the  signal, 
usually  with  a delay  equal  to  an  integer  number  of  sample  spaces.  However, 
delays  of  an  odd  number  of  half-sample  spaces  were  also  investigated,  since 
at  such  values  of  delay  the  extractabillty  of  the  multipath  parameters  from 
the  cepstrum  is  poorest. 

Initially,  it  was  attempted  to  work  with  the  complex  cepstrum,  thus 
leaving  open  the  option  of  eliminating  the  multipath  interference  by  filter- 
ing in  the  cepstrum  domain  ("homomorphic  filtering").  It  was  found,  however, 

as  pointed  out  in  Sec. 2. 3,  that  the  random  fluctuations  in  the  complex  signal 
cepstrum  were  so  great  as  to  mask  the  contribution  of  the  channel.  This  prob- 
lem could  be  overcome  by  working  only  with  the  real  cepstrum.  The  undlstorted 
signal  cannot  be  reconstructed  from  the  real  cepstrum  since  the  phase  spectrum 
has  been  lost;  but  the  multipath  parameters  can  be  estimated  in  the  cepstrum 
(except  for  the  polarity  of  t,  which  must  be  determined  by  other  means)  so 
that  the  adjustment  of  a cancellation  filter  can  be  made  to  follow  the  running 
sequence  of  estimates  obtained  in  this  manner. 

On  the  basis  of  the  data  obtained  with  512-point  sample  sequences,  with 
a ■ 0.2  a multipath  detection  probability  greater  than  0.95  was  estimated 
over  a range  of  t from  about  1 ms  to  about  3.7  ms  (see  Fig.  2.2  of  Cl 3 ).  In 
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the  vicinity  of  t ■ 2 ms.,  the  estimates  of  a were  within  + .07  with  proba- 
bility 0.95  (see  Fig.  28  of  [1]). 

The  steady-state  signal- to-interf erence  ratio  at  the  output  of  an  ideal 
cancellation  filter  was  determined  as  a function  of  the  error  in  the  esti- 
mates of  the  multipath  parameters,  under  the  assumption  of  a time-invariant 
multipath  characteristic.  Values  of  20  db  signal-to-interf erence  ratio 
w^re  obtained  with  high  probability  over  a wide  range  of  a (see  Figs.  27  and 
28  of  [1)). 

The  effect  of  varying  the  length  of  the  original  sequence  used  for  cep- 
strum  processing  was  also  explored  [6].  With  the  type  of  signal  described 
above,  it  was  found  that  the  use  of  128-point  sequences  resulted  in  excessive 
errors,  but  that  useful  rer  .ts  could  be  obtained  by  processing  256-polnt 
sequences . 

All  these  results  were  obtained  under  idealized  conditions  as  follows: 
no  additive  noise  was  considered; 

no  timing  jitter,  sampling  error,  and  other  effects  pertaining  to 
the  analog  waveform  and  A/D  conversion  were  considered; 

only  a single,  discrete,  fixed  multipath  component  was 
considered; 

negligible  round-off  error  was  introduced  into  the  computations 
by  the  635  computer. 
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4.  THE  CURRENT  EFFORT 


4.1  Objectives 

The  earlier  work,  summarized  in  the  preceding  section,  indicated 
that  considerable  reduction  in  multipath  distortion  can  be  achieved  with  the 
aid  of  cepstrum  processing.  However,  these  results  were  obtained  under  very 
idealized  conditions.  The  next  phase  of  work,  reported  herein,  was  therefore 
initiated  with  the  aim  of  developing  and  evaluating  a test  setup  which  would 
correspond  more  closely  to  a practical  signal  processor.  Specifically,  it 
was  desired  to  process  actual  modem  waveforms  (inputs  to  a receiving  modem). 
This  objective  can  be  divided  into  two  stages: 

a)  modem  waveforms  are  locally  generated,  with  multipath  distortion 
artificially  introduced; 

b)  modem  waveforms  are  demodulated  test  transmissions  which  have  been 
sent  over  a long  distance  HF  loop  or  other  type  of  multipath  channel. 

Furthermore,  in  order  to  operate  within  the  typical  limitations  of  a 
practical  special  purpose  digital  processor  which  is  dedicated  to  on-line 
cepstrum  processing  (e.g.,  small  word  size,  sm;  11  memory),  the  necessary 
processing  operations  were  to  be  implemented  on  a small  computer.  The  aim 
was  to  use  the  resulting  test  setup  first  on  signals  incorporating  a single 
discrete  multipath  component  in  order  to  examine  the  extractability  of  multi- 
path  Information  via  the  cepstrum,  for  different  types  of  modems.  This  could 
then  be  extended,  depending  on  available  time,  to  more  complicated  multipath 
structures,  and  to  an  examination  of  the  effects  of  additive  noise. 

Because  the  usefulness  of  the  cepstrum  for  extracting  multipath  informa- 
tion depends  strongly  on  the  nature  of  the  signal,  another  phase  of  the 
current  effort  had  the  aim  of  identifying,  or  designing,  types  of  waveforms 
which  are  best  suited  for  such  processing. 

4.2  Accomplishments 

4*.  2.1  On-Line  Cepstrum  Analysis 

The  development  of  a basic  on-line  cepstrum  processing  capa- 

bilitv  was  completed,  using  a PDP-8I  computer  at  DICEF*.  The  computer 

* 

Digital  Communication  Experimental  Facility,  RADC. 
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programs,  and  the  overall  teat  setup  are  described  in  Appendix  A.  It  turned 
out  that  a major  part  o£  the  overall  effort  vent  into  writing  and  debugging 
the  PDP-8  programs.  This  was  mainly  because  all  programming  was  done  in 
machine  language,  in  order  to  minimize  core  storage  requirements  and  thus  to 
be  able  to  accommodate  all  programs,  and  also  provide  adequate  space  for  data, 
within  the  4K  core. 

For  these  reasons,  the  debugging  of  the  various  routines  was  completed 
only  near  the  end  of  the  project  period,  so  that  there  remained  not  much  time 
to  collect  test  data. 

Problems  with  equipment  also  caused  some  delay.  One  of  these  was  a con- 
trol circuit  failure  in  an  incremental  read-write  magnetic  tape  unit  which 
had  been  connected  to  the  PDP-8I.  Some  time  had  been  spent  to  program  the 
PDP-8I  to  work  with  this  tape  deck,  in  order  that  programs  could  be  loaded 
from  tape,  and  data  recorded  on  the  tape.  The  tape  deck  had  also  been  viewed 
as  an  augmentation  of  available  storage  capacity.  Thus,  additional  programs 
(for  instance,  detection  of  two  multipath  components)  could  have  been  called 
into  core  when  needed,  and  therefore  would  not  require  a permanent  core  allo- 
cation. However,  the  failure  which  occurred  could  not  be  corrected  in  time 
for  the  tape  deck  to  be  used  on  this  project. 

4.2.2  Test  Results 

The  completed  on-line  cepstrum  analysis  and  multipath  param- 
eter estimation  system  (as  described  In  Appendix  A)  was  operated  with  signals 
derived  from  three  different  data  modems,  as  described  in  Appendix  C.  A test 
consisting  of  an  extensive  series  of  detection  runs  was  made  only  with  the 
Codex  9600  modem.  However,  the  limited  data  which  was  obtained  with  the 
other  tvo  modems  .indicates  that  the  Codex  results  can  be  taken  as  representa- 
tive for  the  performance  of  all  three  types  of  modems.  The  data  from  this 
test  is  presented  in  Section  C.2.1  and  is  analyzed  in  Appendix  D.  The  test 
was  conducted  with  the  following  multipath  parameters: 
o ■ 0.5,  t ■ 2.9  ms. 

The  sampling  rate  was  approximately  7000  samples/sec.,  which  made  tne  multi- 
path  delay  approximately  equal  to  20.3  sample  spaclngs. 

The  results  of  the  Codex  detection  test  were  disappointing  for  two 


reasons  which  will  be  discussed  in  turn,  namexy: 

a)  the  wide  spread  in  the  estimate  of  n;  and 

b)  the  large  number  of  detection  errors. 

The  standard  deviation  for  the  set  of  a-estimates  (using  only  those  esti- 
mates which  resulted  from  correct  detections)  is  0.102;  furthermore,  this 
value  would  be  even  larger  (by  about  10Z)  if  a correction  for  multipath  cep- 
strum  depression  had  been  applied  to  the  data.  In  order  to  be  able  to  esti- 
mate the  contribution,  to  the  variability  of  the  a-estlmates,  of  various 
errors  Introduced  in  processing,  computations  of  these  errors  are  carried  out 
in  Appendix  D (based  on  the  analysis  of  errors  in  Appendix  B).  It  is  shown 
there  that  the  largest  source  of  error  is  the  timing  jitter  which  was  inad- 
vertently introduced  into  the  sampling  process.  Nevertheless,  based  on  the 
approximations  of  Appendix  B and  Appendix  D,  a rather  small  estimate  of  the 
total  processing  error  was  obtained.  As  is  pointed  out  in  Sec.  D.2,  however, 
this  estimate  is  very  likely  to  be  on  the  low  side. 

In  summary,  the  data  and  analysis  suggests  that,  although  errors  intro- 
duced in  processing  have  produced  a noticeable  effect  in  the  cepstrum,  this 
effect  does  not  account  for  the  major  part  of  the  variability  in  a-estlmates. 
It  must  therefore  be  concluded  that  >t  good  part  of  the  variability  is  inher- 
ent in  the  particular  types  of  modem  signals  which  were  used. 

Turning  now  to  the  large  number  of  detection  errors  (25  out  of  a total 
of  201  multipath  estimates),  it  is  clear  that  the  above  comments  apply  here 
also.  Both  of  the  plots  of  signal  cepstrum  statistics  (Figs.  Cl  and  C7) 
imply  a rather  large  variability  „Z  cne  signal  cepstrum.  The  auestlon  as 
to  the  extent  to  which  this  is  du':  to  errors  introduced  in  processing  is 
•"gain  addressed  by  Appendix  D.  On  the  other  hand,  using  the  standard  devia- 
tion of  the  (noisy)  signal  cepstrum,  a computation  discussed  in  Section  D.l 
yields  a much  smaller  probability  of  wrong  detection  than  actually  obtained. 
Reasons  for  this  discrepancy  are  listed  in  that  Section.  The  question  arises 
whether  processing  errors  significantly  affect  the  shape  of  the  cepstrum  dis- 
tribution and  in  this  way  contribute  disproportionately  to  the  occurrence  of 
detection  errors.  This  has  not  been  investigated. 
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4.2.3  Analysis  of  Other  Factors  Affecting  Practical  Application 


The  tests  described  In  the  preceding  Section  were  of  limited 
scope  in  a number  of  ways: 

a)  Only  a single  multipath  component  was  used  of  fixed  magnitude  and 
delay. 

b)  The  overall  test  was  carried  out  strictly  at  baseband. 

c)  No  actual  multipath  cancellation  was  performed. 

The  relatively  simple  algorithm,  described  in  Sec.  2.4,  for  extracting 
the  multipath  parameters  can  no  longer  be  used  if  the  presence  of  more  than 
one  multipath  component  is  anticipated.  An  algorithm  applicable  to  the  use 
of  two  multipath  components,  that  is,  a channel  with  the  impulse  response 

h(t)  - S(t)  + a$(t  - r^)  + 66 (t  - t 2),  (4-1) 

is  presented  in  Appendix  E.  From  Section  E.l  it  is  clear  that  multipath  de- 
tection, and  parameter  estimation,  under  the  assumption  that  two  multipath 
components  are  present,  is  much  more  involved  than  in  the  case  where  only  a 
single  multipath  component  is  assumed.  Not  only  must  several  cepstrum  peaks 
be  taken  into  account  (whose  magnitudes  are  not  necessarily  linearly  related 
to  the  multipath  coefficients),  but  rather  fine  discrimination  of  their  rela- 
tive magnitudes  may  be  required,  as  is  evidenced  by  Fig.  El.  This  means  that 
the  results  obtained  with  this  algorithm  are  much  more  affected  by  signal 
cepstrum  variability  and  additive  noise  than  those  obtained  in  the  case  of  a 
single  multipath  component.  Extraction  of  multipath  parameters  for  even  more 
complicated  multipath  structures  appears  to  be  out  of  the  question,  except  in 
special  cases.  In  such  cases,  multipath  cancellation  by  homomorphic  filtering 
may  be  possible,  provided  that  an  effective  separation  of  signal  cepstrum  and 
multipath  cepstrum  can  be  arranged. 

In  appendix  G,  the  analvsis  was  extended  to  HF  signals,  using  SSB  and 
DSB  modulation.  It  is  shown  there  (Eq.  G-9)  that  ordinary  SSB  demodulation 
of  a multipath  signal  results  in  a received  baseband  signal  in  which  the 
multipath  appears  in  a more  complicated  way,  involving  the  Hilbert  transform. 
Complex  demodulation,  on  the  other  hand,  gives  the  simple  baseband  signal 
(G-13)  which,  however,  is  complex.  Due  to  the  fact  that  the  received  multi- 
path  amplitude  coefficient  a'  is  then  also  complex,  cepstrum  analysis  will 
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have  to  be  applied  to  the  complex  signal  rather  than  only  Its  real  part.  An- 
other factor  to  be  considered  Is  the  fact  that  the  argument  of  the  coefficient 
a'  Is  -u)qt  (Eq.  G-14) , where  uQ  Is  the  carrier  angular  frequency,  so  that  It  is 
very  sensitive  to  small  changes  In  t.  These  observations  indicate  that  some 
modifications  of  the.  present  processing  scheme  are  required  in  order  to  make 
it  applicable  to  actual  HF  signals,  even  if  only  a single  multipath  component 
is  assumed. 

In  order  to  assess  the  overall  performance  of  a multipath  processor,  the 
performance  of  the  actual  cancellation  scheme  used  must  be  taken  into  account. 
In  [1],  signal- to- interference  ratio  at  the  canceller  output  was  plotted  in 
Fig.  27  as  function  of  a and  a,  for  the  case  of  an  ideal  recursive  cancella- 
tion filter  operating  in  steady  state,  and  without  discrepancy  between T and 
r.  In  practice,  canceller  operation  deviates  from  these  conditions,  particu- 
larly because  starting  transients  will  be  associated  with  each  new  set  of 
parameter  values  fed  to  the  canceller.  Also,  stability  problems  arise  with  a 
recursive  filter  which  may  make  its  use  undesirable.  An  alternative  possibil- 
ity is  to  carry  out  multipath  cancellation  in  the  DFT  domain.  Analogous  to 
the  starting  transient  of  a recursive  filter,  this  results  in  a residual  error 
which  is  greatest  near  the  beginning  of  the  reconstructed  waveform.  This 
erfect  can  be  compensated  for  by  dropping  points  off  the  leading  end  of  the 
reconstructed  waveform.  As  an  example,  using  N-256  samples  of  a modem  signal 
with  multipath  at  a **  0.5,  t * 20  T/N,  and  basing  the  cancellation  on  the  ex- 
act multipath  parameters,  an  r.ra.s.  signal-to-interf erence  ratio  over  all  256 
points  of  only  14  db  was  obtained.  However,  a 1 db  improvement  accompanied 
every  4 points  dropped,  with  S/I  * 30  db  when  only  196  points  were  retained. 

4.2.4  Signal  Design  [ 3 ] , [8 ] 

In  addition  to  the  computer  programming  and  modem  tests 
which  have  been  discussed  in  the  preceding  Sections,  another  major  area  of 
activity  has  been  that  of  signal  design.  Whereas  the  tests  with  the  various 
data  modems  were  carried  out  with  the  idea  that  multipath  compensation  would 
be  performed  strictly  as  supplementary  to  the  normal  operation  of  an  existing 
system,  the  question  addressed  in  the  signal  design  effort  was  this:  What  is 

the  best  (or  a good)  transmitted  waveform  to  be  used  if  a complete  system 
design  is  to  be  pursued.  In  the  following  paragraphs  is  outlined  the  general 


approach  which  has  been  taken  to  the  problem  of  signal  design  for  multipath 
channels. 


It  is  assumed  that  the  received  signal  is  disturbed  by  discrete  multi- 
path  components  which  can  be  assumed  constant  over  one  observation  period  (a 
sequence  of  N samples),  and  that  the  multipath  parameters  needed  for  cancel- 
ling or  reducing  the  multipath  distortion  are  estimated  either  in  part  or  in 
total  by  cepstrum  calculations. 

The  signal  design  problem  is  therefore  to  minimize  cepstrum  variability 
due  to  the  transmitted  signal,  in  order  that  the  multipath  parameters  become 
as  clearly  discernible  in  the  cepstrum  domain  as  possible.  For  instance, 
the  signal  cepstrum  might  be  made  zero,  or  close  to  zero,  over  that  portion 
of  the  t'-domaln  which  contains  the  detection  window  - the  portion  needed  for 
detecting  multipath  and  estimating  the  multipath  parameters. 


A signal  with  a small  real  cepstrum  is  one  whose  Fourier  transform  has 
approximately  constant  magnitude  over  the  frequency  range  of  interest.  In 
order  to  transmit  information  by  means  of  such  signals,  the  transmitted  in- 
formation must  therefore  reside  essentially  in  the  phase  spectrum  of  the 
signal.  This  can  be  accomplished  by  generating,  within  time  intervals  of 
length  T,  transmitted  signals  x(t)  as  a convolution. 


x ( t)  - xp(t)*xn(t) , (4-2) 

where  x (t)  is  a pulse  shape  giving  rise  to  a suitably  flat  spectrum,  and  x 
p m 

is  an  all-pass  modulating  sequence,  x (t)  ■ E a,5(t  - IT  ) . This  modulating 

m l m 

sequence  is  obtained  through  a suitable  transformation  of  the  data  sequence 
to  be  transmitted.  Such  a transformation  is  obtained  if  the  message  sequence 
is  used  directly  to  define  a nhase  sesuenc"  in  the  frequency  domain.  Thus, 
let 


dx  (n)  » 1,  n - 0,  ...,  N-l 

d^n)  - ej2*Q(n),  > (4"3> 

m 

where  Q(n)  represents  a data  sequence.  Then,  after  forming  dx(n)  ■ 

dx  (n)*dx  (n) , the  inverse  DFT  can  be  applied  and  a time  waveform  produced  by 

pm 

a suitable  smoothing. 
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The  details  require  a careful  examination  of  suitable  ways  of  comput- 
ing the  eepstrum,  and  of  various  properties  associated  with  different  ver- 
sions of  cepstra.  For  this  reason,  various  definitions  of  real  and  complex 
cepstra  were  established,  compared,  and  studied  in  detail,  as  presented  in 
[3].  This  Included  modifications  and  extensions  of  previous  developments 
to  Include  practical  considerations  such  as  bandwidth  and  time-duration  lim- 
itations. It  could  then  be  shown  that  the  proposed  signaling  scheme  does  not 
require  sampling  synchronization  between  the  receiver  and  transmitter.  How- 
ever, frame  synchronization  is  required;  that  is,  a received  data  sequence 
must  be  derived  from  only  a single  transmitted  waveform  of  the  type  (4-2) . 

A careful  examination  of  this  synchronization  requirment,  and  of  the  effects 
of  overlap  between  successive  transmissions  has  not  been  made  as  yet. 

4.2.5  Effects  of  Additive  Noise 

The  effect  in  the  eepstrum  due  to  additive  noise  in  the  time  domain 
was  investigated  in  a related  project  [2].  This  work  was  carried  out  by 
means  of  computer  simulation,  using  APL.  A brief  summary  is  presented  in 
Appendix  F.  The  empirical  formula  (F-2)  presented  there  was  applied  to  the 
error  analysis  in  Appendix  D. 

No  tests  with  additive  noise  were  performed  with  the  experimental 
setup  and  PDP-8  eepstrum  processor. 


5.  CONCLUSIONS 

a)  The  basic  objective  of  implementing  an  on-line  cepstrum  pro- 
cessor, for  detecting  the  presence  of  a single  multipath  component  in  a 
baseband  signal  and  estimating  its  magnitude  and  delay,  was  accomplished. 

b)  Only  a very  limited  amount  of  testing  was  possible.  However, 
since  the  test-setup  and  computer  programs  continue  to  be  available,  it 
would  be  possible  to  run  an  extended  series  of  tests  under  a variety  of 
controlled  conditions,  and  thus  obtain  more  comprehensive  and  reliable  per- 
formance data. 

c)  On  the  basis  of  the  limited  test  data  which  was  obtained,  a 
rather  marginal  performance  was  observed.  A major  source  of  error  (timing 
jitter  in  sampling)  was  not  recognized  until  after  the  project  terminated. 
With  additional  time,  its  correction,  as  well  as  a number  of  other  adjust- 
ments (such  as  data  window  shaping,  or  varying  the  domain  of  the  logarithm 
function)  would  have  produced  some  improvement  in  performance.  However, 
based  on  the  analysis  of  the  effects  of  various  sources  of  error,  such 
modifications  cannot  be  expected  to  result  in  a large  reduction  in  the 
variability  of  multipath  parameter  estimates,  below  that  observed  in  the 
tes*  data  of  Appendix  C.  In  other  words,  the  modem  signals  themselves 
produce  considerable  variability  in  the  cepstrum. 

d)  When  this  type  of  cepstrum  processor  is  operated  as  part  of  a 
practical  communication  link,  error  in  multipath  cancellation  will  arise 
not  only  from  the  variability  of  the  signal  cepstrum,  but  from  additive 
noise,  from  fluctuations  of  the  multipath  parameters  with  time,  and  from 
transients  inherent  in  the  cancellation  system.  This  suggests  that  useful 
operation  will  hardly  be  feasible  unless  the  signal  cepstrum  variability 
is  carefully  controlled,  which  is  achievable  through  proper  signal  design. 

e)  The  requirement  of  low  signal  cepstrum  variability  becomes  even 
more  urgent  when  the  results  presented  in  Appendix  E are  taken  into  account. 
A practical  HF  system  would  have  to  be  able  to  compensate  for  at  least  two 
multipath  components;  and  in  that  ease,  parameter  estimation  requires  a 
rather  fine  discrimination  between  cepstrum  peaks.  Use  of  a combination 

of  techniques  (e.g.,  cepstrum  analysis  plus  autocorrelation  analysis)  may 
be  a viable  strategy  for  achieving  satisfactory  performance. 
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APPENDIX  A.  DESCRIPTION  OF  PDP-8  PROGRAMS  AND  TEST  SETUP 

A.l  General  Description  of  Programs 

The  basic  functions  which  have  been  programmed  and  are  available  within 
the  4K(4096-word)  core  memory  of  the  PDP-8  are  itemized  below.  Fig.  A1  shows 
the  core  allocation  schedule  by  memory  page.  (Page  numbers  are  base-8) . 

a)  Storage  of  a 512-point  sequence  of  data  samples  (x^(k) • 

(Memory  pages  7-12;  see  also  Section  A. 2.1.) 

b)  Application  of  a data  window,  and  overlapping  of  the  512-point  se- 

255 

quence  into  a 256-point  sequence.  That  is,  a new  sequence  (x(k) 
is  generated,  where  x(k)  * ■|'[w(k)x^(k)  + w(k  + 256)x^(k  + 256)]. 
(Memory  page  13;  see  also  Section  A. 2. 2.) 

c)  Computation  of  the  DFT  of  (x(t)},  giving  {d^(n)}.  For  this  prrpose, 
the  complex  FFT  routine  from  the  DECUS  Program  Librar  * is  used.  Al- 
though the  cepstrum  computation  as  it  is  now  carried  out  requires  only 
the  real  FFT  (discrete  Fourier  transform  of  a real  sequence),  the 
complex  FFT  was  used  in  developing  the  program  package  in  order  to 
provide  the  option  of  carrying  along  the  imaginary  part.  (Memory 
pages  2-6;  the  result  of  this  computation  appears  in  pages  7-10  (real 
part)  and  pages  17-20  (*maginary  part).  See  also  Section  A. 2. 3.) 

d)  Computation  of  the  logarithm  of  the  magnitude  of  (dx(n)}. 

(Memory  pages  30-31;  see  also  Section  A. 2.4.) 

e)  Computation  of  the  inverse  DFT  (using  again  the  complex  FFT  routine) 
to  yield  the  real  cepstrum  (cx(k)}.  (See  also  Section  A. 2. 5.) 

f)  Computation  of  the  mean,  maximum  and  minimum  values  at  each  of  the 
first  40  cepstrum  points,  for  16  successively  computed  cepstra. 

(Memory  pages  23-26;  see  also  Section  A. 2. 7.) 

g)  Threshold  detection  of  a single  multipath  spike  in  the  cepstrum,  with 
a printout  of  estimated  amplitude  and  position  of  the  spike. 

(Memory  pages  14-16  and  22;  see  also  Section  A. 2. 6.) 

The  main  calling  program  for  carrying  out  the  cepstrum  computation  is 
stored  in  Page  1.  This  includes  a routine  for  acquiring  a sequence  of  signal 
samples  via  A/D  conversion.  There  is  also  a provision  for  oscilloscope  dis- 
play of  stored  data,  via  D/A  conversion. 
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Page  0 is  reserved,  in  the  usual  manner,  for  constants  and  pointer  ad- 
dresses. Pages  34-37  contain  the  loaders  and  an  Edit  Routine  devised  by 
Capt.  Russell  Lemon  of  RADC.  The  Edit  Routine  facilitates  the  editing  of 
programs  and  provides  for  various  other  manipulations.  Finally,  pages  32-33 
contain  routines  for  reading  and  writing  onto  a magnetic  tape  unit  which  is 
one  of  the  I/O  devices  incorporated  in  the  test  setup  at  DICEF. 

A. 2 Detailed  Description  of  Processing  Steps. 

A. 2.1  Data  Acquisition 

Assume  that  a signal  waveform  has  been  connected  to  the  A/D  converter 
input,  with  peak  values  not  exceeding  the  dynamic  range  of  the  converter. 
Under  program  control,  A/D  conversion  is  Initiated  and  continued  until  a se- 
quence of  512  signal  samples  (at  a preselected  sampling  rate)  has  been  stored 
in  memory  locations  1600  through  2577  (pages  7-12) . Data  acquisition  occurs 
under  control  of  the  Main  Calling  Program. 

A. 2.2  Data  Preparation 

Before  Fourier  transformation,  a data  window  is  applied,  and  the  512- 
point  sequence  is  split  and  overlapped  into  a single  256-point  sequence.  The 
data  window  which  is  used  results  in  a raised-cosine  weighting  of  the  first 
and  last  25  points  of  the  data  sequence,  similarly  to  the  window  used  in  [1]. 
Specifically,  if  (x^k)}  denotes  the  512-point  sample  sequence  (k  - 0,1  ..., 
511),  then  {x^(k)}gets  multiplied  by  the  sequence  (w(k)},  where 

/ 0.5(1  - cos  it),  k - 0,  1 25 

w(k)  - < 1,  k - 26 486  (A-l) 

[ 0.5(1  - cos  '^5i28~  k - 487,  488 511. 

(The  particular  form  for  the  argument  of  the  cosine  resulted  from  the  fact 
that  the  trigonometric  table  of  the  FFT  routine  is  used  to  obtain  the  cosine 
values.) 

Splitting  and  overlapping  the  weighted  512-point  data  sequence  results 
in  a new  sequence  (x(k)}  (k  - 0,  1,  ...,  255),  where 

x(k)  - [w(k)xx(k)  + w(k  + 256)xx(k  + 25b)].  (A-2) 
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The  factor  is  needed  to  avoid  overflow.  These  operations  are  carried  out 
by  the  subroutine  WINDOW,  which  deposits  the  final  result  {x(k)}  in  memory 
locations  1600  through  2177  (pages  7 and  10) . 

A.2.3  Discrete  Fourier  Transform 

The  discrete  Fourier  transform  of  the  sequence  (x(k)}  is  computed  using 
the  FFT  routine  for  complex  data  which  is  available  through  the  DECUS  Program 
Library  [5].  The  result  is  a complex  sequence  {dx(n)}  (n  ■ 0,  1,  ...»  255). 

It  Includes  a scale  factor  (a  power  of  2)  which  is  automatically  assigned  by 
the  FFT  routine  in  order  to  utilize  the  limited  word  rize  as  effectively  as 
possible,  thus  giving  a quasi  floating-point  effect.  However,  this  scale 
factor  becomes  merely  an  additive  constant  upon  taking  the  logarithm,  and 
therefore  contributes  only  to  the  first  cepstrum  point,  c (0).  Since  c (0) 

X X 

is  * no  significance  in  later  processing  and  in  fact  will  be  arbitrarily  re- 
assigned later,  it  is  not  necessary  to  carry  along  the  scale  factor  in  question. 

The  FFT  routine  operates  on  data  stored  in  memory  pages  7,  10,  17  and  20. 

It  considers  pages  7 and  10  the  real  part  of  the  input  data,  and  pages  17  and 
20  as  the  Imaginary  part.  However,  in  the  present  program  package  only  real 
inputs  are  provided  for,  and  only  the  real  cepstrum  is  computed,  so  that  the 
FFT  routine  is  used  only  for  real  data.  Therefore,  pages  17  and  20  are  set 
to  zero  prior  to  execution  of  the  FFT  routine.  The  real  part  of  the  result 
of  the  FFT  computation  resides  in  pages  7 and  10,  the  imaginary  part  in  pages 
17  and  20. 

A.2.4  Logarithm  of  Absolute  Magnitude 

In  the  LOG  subroutine,  first  the  squared  magnitude  of  {d^Cn)}  is  computed, 
by  squaring,  and  then  summing,  the  real  and  imaginary  parts  for  each  n.  (The 
square  root  is  not  taken,  since  it  amounts  to  merely  a scalar  factor  after 
taking  the  logarithm,  and  this  becomes  part  of  the  scaling  performed  in  the 
cepstrum  domain.)  The  additional  12  digits  which  appear  in  the  MQ  (multi- 
plier quotient)  register  after  squaring  are  carried  along  and  used  in  the 
\ 

logarithm  computation,  so  that  no  truncation  error  is  incurred  in  this  step. 

The  logarithm  is  computed  using  a look-up  table  for  the  seven  most  sig- 
nificant digits:,  which  yields  the  mantissa,  and  the  number  of  leading  zeros 


24 


/ 


i 


gives  the  characteristic.  This  is  Illustrated  in  Fig.  A2.  It  is  the  binary 
logarithm  which  is  actually  computed,  so  that  another  scale  factor  must  be 
applied  in  the  cepstrum  domain  in  order  to  achieve  the  same  result  as  would 
be  obtained  with  the  natural  logarithm. 

The  result  of  the  logarithm  computation,  the  sequence  (i^Cn)},  is  stored 
in  pages  7 and  10. 
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Fig.  A2.  Logarithm  Computation 
A. 2.5  The  Real  Cepstrum 

The  real  cepstrum  is  obtained  by  performing  the  inverse  DFT  on  the  real 
logarithm,  using  the  FFT  routine  in  the  manner  described  in  Section  A. 2. 3. 

In  order  to  achieve  maximum  resolution  in  the  result  of  the  FFT  computa- 
tion, the  output  sequence  should  not  have  excessively  large  peaks  so  that 
advantage  can  be  derived  from  the  automatic  scaling  introduced  by  the  FFT 
routine  (see  Sec.  A. 2. 3).  Usually,  in  cepstrum  computation,  a very  large 
value  of  c (0)  is  obtained,  due  to  the  d.c.  component  in  {l  (n)}.  In  order 
to  avoid  this  large  peak  in  the  cepstrum,  the  log  spectrum  has  to  be  shifted 
to  make  its  d.c.  component  zero,  before  the  inverse  DFT  is  computed.  This 
step  is  carried  out  in  the  LOG  routine. 

The  scale  factor  introduced  in  the  FFT  routine  during  the  inverse  DFT 
computation  must  be  corrected  for  in  the  cepstrum. 


The  processing  steps  described  in  Sections  A. 2.1  through  A. 2. 5 are  all 
executed  under  control  of  the  Main  Calling  Program. 


A. 2.6  Detection  of  Multipath 


The  purpose  of  the  Detection  Routine  is  to  decide  whether  or  not  multi- 
path  is  present  and,  if  the  decision  is  that  multipath  is  present,  to  estimate 
the  multipath  amplitude  and  delay  (Assuming  a single  discrete  multipath  com- 
ponent). The  Detection  Routine  can  be  executed  Independently,  or  it  can  be 
brought  under  the  control  of  the  Main  Calling  Program. 

The  Detection  Routine  consists  of  the  Detection  Calling  Program  and  the 
following  subroutines: 

a)  Threshold  Comparison  Routine 

b)  a-Estimation  Routine. 

The  Detection  Calling  Program  applies  all  the  required  scale  factors  to 
the  real  cepstrum  sequence  residing  in  pages  7 and  10,  except  for  the  factor 
In  2 ■ .693  which  converts  the  binary  logarithm  to  the  natural  logarithm  (see 
Sec.  A. 2. A).  Let  the  result  obtained  after  scaling  be  denoted  ( k * 0, 

1,  ...,  255).  The  program  then  subtracts  from  the  first  100g  ■ 64  cepstrum 

points  the  sequence  of  mean  values  {c^}  (i  - 0,  1 63)  which  are  stored 

in  locations  5000  - 5077  (See  Sec.  A. 2. 7)  and  deposits  the  resulting  sequence 
(c’}  * {c^}  - {c^}  (i  - 0,  1,  ...,  63)  in  locations  1600  - 1677.  The  Thresh- 
old Comparison  Routine  is  then  executed. 

A. 2.6.1  Threshold  Comparison 

The  sequence  {c^}  is  compared  against  the  two  threshold  sequences  {6*} 
and  {0^}  (see  Sec.  A. 2. 7)  which  are  stored  in  locations  5100  - 5277.  This 
results  in  the  sequence  {q^}  (i  » 0,  1,  ...,  63)  defined  as  follows: 

q.  ■ 0 for  each  i such  that  0 . <,  c'  £ 0* 

q^  ■ cj^  - 0^  for  each  i such  that  c^  > 0^ 

q^  ■ 0^  - c|  for  each  i such  that  < 8^,  . 

{q^}  is  entered  into  locations  5600  - 5677.  The  largest  number  in  the  se- 
quence {q^},  max  q^,  is  then  determined,  and  the  corresponding  subscript  i - i. 
If  max  q^  » 0,  then  the  message  "NO  MP"  (no  multipath)  is  typed  out  on  the 
typewriter  and  execution  terminates. 
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If  max  q^  > 0,  execution  proceeds  on  the  assumption  that  a multipath 
spike  has  been  detected.  On  the  typewriter  the  cepstrum  value  cl  and  its 
storage  location  (1600  + i)  is  typed  out  for  reference.  The  two  'numbers  cj 
(j  ■ i+l,*i-l)  are  then  examined  to  determine  whether  either  one  satisfies 
the  following  three  requirements: 

a)  same  polarity  as  cl  , 

b)  q^  > 0 , 

c)  q.  - max  q , (1  - £+1,1-1). 

J l 

If  this  test  is  satisfied  by  either  j-value,  denote  it  j.  The  numbers  j and 
cl  are  stored.  The  program  then  jumps  to  the  a-Estimation  Routine. 

A. 2. 6. 2 a-Estimation  Routine 


If  no  j is  stored,  a is  estimated  from  cl  only.  That  is,  estimation  of  a 

1 . 

proceeds  on  the  assumption  that  t Is  exactly  equal  to  i sampling  intervals. 


The  estimate  of  a is  then 

t 


.693  c 

n7 T 


1 


(A-3) 


and  the  estimate  of  t is  t * i/f  , where  f **  T/N  is  the  sampling  frequency. 

» s s 


On  the  other  hand,  if  j has  been  stored,  then  the  a-estimate  is  based  on 
boph  cl  and  cl,  and  an  interpolation  between  cl  and  cl  is  performed  as  des- 
cribed in  Sec?  2.4.  That  is,  let 


♦ * lcjl  / i !cj I + |c:|). 


Then  the  a-estimate  is 
Cj  x 2 x .693 

- N sine  $ 
and  the  estimate  of  t is 

T . i±  d - m 

1 f 

8 


(A-4) 


(A-5) 


(The  present  version  of  the  program  does  not  provide  for  explicit  computation 
Of  T.) 


A A 

A typewrite**  printout  is  then  generated,  which  lists  a,  i,  j (i£  it 

exists),  ! cl  I , I cl  I , and  fell  /sine  <$> . Thereafter  the  Detection  Routine  ter- 

i J i 

minates. 

A. 2. 7 Cepstrum  Statistics 

The  Detection  Routine  requires  mean  cepstra  and  thresholds  for  any  given 
type  of  modem  signal.  Using  a separate  Calling  Program,  statistics  can  be 
compiled  for  the  first  100g  ■ 64  cepstrum  points  for  a set  of  16  signal  se- 
quences. This  Calling  Program  cycles  through  the  Main  Calling  Program  16 
times  in  succession,  and  processes  the  first  64  points  of  each  of  the  16 
cepstra  thus  obtained  to  yield 

a)  a mean  cepstrum  (estimate)  for  the  given  signal  source,  fc^}, 
i m 0,  1,  . . . , 63 

b)  the  sequence  of  maximum  values 

c)  the  sequence  of  minimum  values  (negative  maxima) 

The  sequence  {c^}  is  entered  into  locations  5000  - 5077,  the  maxima  in  loca- 
tions 5100  - 5177,  the  minima  in  locations  5200  - 5277.  The  mean  values  are 
used  directly  by  the  Detection  Routine  (Sec.  A. 2. 6).  The  minima  and  maxima 
can  be  used  to  establish  thresholds.  However,  this  has  not  been  programmed; 
Instead,  thresholds  have  been  entered  manually  for  the  various  test  runs. 

A. 3 The  General  Test  Setup 

In  order  to  perform  on-line  cepstrum  analysis  of  communication  signals, 
the  PDP-8I  at  DICEF  was  operated  in  connection  with  other  equipment  as 
illustrated  in  Fig.  A3.  A data  modem  is  driven  from  a transmission  test  set 
which  generates  a long  pseudo-random  binary  modulating  sequence.  The  base- 
band output  of  the  modem  is  returned  to  the  baseband  input  for  demodulation 
and  the  resulting  data  stream  returned  to  the  transmission  test  set,  for  an 
error  check.  This  connection  allows  verification  that  a proper  baseband  out- 
put is  derived  from  the  modem. 

The  modem  output  is  also  fed  through  the  DICEF  Hireline  Simulator  which 
is  programmed  to  introduce  the  desired  multipath  interference.  (For  cepstrum 
computations  of  the  undistorted  modem  signal,  the  wireline  simulator  is  by- 
passed). Since  the  conversion  range  of  the  A/D  converter  is  + 10  V,  the 
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Wireline  Simulator  output  is  passed  through  variable  attenuation  and  an 
amplifier. 

A/D  conversion  takes  place  under  program  control.  That  is,  conversion 
of  a sequence  of  samples  is  initiated  when  the  computer  is  ready  to  accept  the 
data.  A square  wave  generator  connected  to  the  A/D  converter  determines  the 
sampling  rate. 

Oscilloscope  observation  of  stored  data,  such  as  a cepstrum  sequence,  is 
possible  through  a D/A  routine  w?iich  reads  out  the  desired  data  sequence  re- 
spectively to  a D/A  converter.  The  oscilloscope  is  connected  to  the  D/A  con- 
verter output. 


Fig.  A3  Test  Setup 
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APPENDIX  B.  SOURCES  OF  ERROR 


B. 1 Accumulated  effect  of  truncation  error  In  the  cecstrum 

The  accumulated  effect  of  truncation  errors  which  occur  in  the  FFT  and 
logarithm  computations  was  examined  by  computing  the  cepstrum  of  a signal 
sequence  both  via  PDP-8  and  also  via  APL\370.  For  this  purpose,  a 512  point 
sample  sequence  derived  from  a USC-9  modem,  with  multipath,  was  first  sub- 
jected to  the  WINDOW  routine  (that  is,  the  data  window  was  applied  and  the 
sequence  was  folded  over  into  256  points.)  The  resulting  256-point  sequence, 
in  4-digit  octal  representation,  was  typed  into  the  APL  terminal,  converted 
to  decimal  form,  and  its  cepstrum  computed.  The  cepstrum  was  also  computed 
via  PDP-8.  The  APL  computation  yields  the  cepstrum  with  negligible  error 
due  to  truncation.  Comparison  of  the  two  cepstra  thus  obtained  was  accom- 
plished by  again  copying  the  PDP-8  result  into  APL,  and  computing  the  differ- 
ence of  the  two  cepstra.  The  resulting  error  sequence  (excluding  the  initial 
point)  was  characterized  by  the  following  parameter  values: 

Peak  error:  0.755 

RMS  error:  0.267 

These  values  must  be  divided  by  N/2  - 128  in  order  to  relate  them  to  a,  the 
estimated  multipath  amplitude.  Thus,  if  the  peak  error  of  0.755  were  to 
occur  At  exactly  that  point  In  the  cepstrum  where  the  largest  multipath  spike 
is  located,  and  assuming  that  the  multipath  delay  t is  an  exact  multiple  of 
the  sample  spacing,  then  this  peak  truncation  error  will  contribute  an  error 
of  +0.755/128  a +0.006  to  o.  The  r.m.s.  error  in  jj  due  to  truncation,  on  the 
basis  of  this  single  test,  is,  therefore,  0.267/128  * 0.002,  which  is  small. 

This  result  must  be  accepted  with  some  reservations,  however.  The 
amplitude  distribution  of  (d^(n)},  the  DFT  sequence,  can  vary  noticeably 
for  different  input  sequences  from  the  same  source.  This  distribution 
affects  the  magnitude  of  the  truncation  error,  that  is,  more  truncation  error 
is  associated  with  the  very  small  values  of  (dy(n)}. 

B.2  Error  introduced  in  Data  Preparation 

The  truncation  error  check  described  in  the  preceding  section  does  not 
account  for  error  introduced  during  data  preparation  using  the  WINDOW  routine 
(sec  Sec.  A. 2. 2).  Truncation  errors  are  introduced  in  the  representation  of 
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cosine  values,  the  multiplication,  and  the  scaling.  However,  the  errors 
due  to  the  truncated  representation  of  the  cosine  values  are  not  variable, 
since  a fixed  look-up  table  is  used.  Therefore,  they  do  not  contribute  to 
the  variability  of  the  u-estimate.  The  same  statement  applies  to  the  w(k)'e 
in  (A-l). 

Scaling  and  overlapping, as  indicated  in  (A-2) , results  in  an  error  dis- 
tribution at  each  sample  point  as  shown  in  Fig.  Bl,  which  has  a variance  of 
1/8.  At  each  of  the  first  and  last  25  points  cf  the  resulting  256-point 
sequence,  one  of  the  two  numbers  being  added  is  the  result  of  a multiplication 
and,  therefore,  has  sustained  an  additional  truncation  error.  The  error  dis- 
tribution at  these  points  is  therefore  approximately  of  the  form  shown  in 
Fig.  B2,  and  has  a variance  of  7/48.  But  this  value  only  applies  at  approxi- 
mately 20%  of  the  256  points  and  does  not  differ  greatly  from  1/8,  so  that 
the  total  error  introduced  in  conditioning  the  input  sequence  has  a variance 
of  approximately  1/8.  The  corresponding  error  standard  deviation  is 
VTTs  - 0.35. 

This  error  represents  additive  noise  at  the  input  to  the  cepstrum 
processor.  Its  effect  on  the  cepstrum  can  be  estimated  by  the  method  of 
Appendix  F. 


1 


.0  1 

Figure  Bl  Figure  B2 

Error  Distribution  Error  Distribution 

B . 3 Effect  of  Signal  Quantization 

The  effect  of  signal  quantization  on  the  cepstrum  was  explored  by  com- 
puting the  cepstrum  of  the  same  signal  sequence  several  times,  but  each  time 
dropping  an  additional  digit  in  the  binary  representation  of  each  input  signal 
sample.  Thus,  a 256-point  signal  sequence  was  recorded  on  paper  tape,  and  its 
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cepstrum  computed.  The  same  signal  sequence  was  then  again  read  into  the 
PDP-8,  a right-shift  of  one  digit  applied,  and  the  cepstum  computed.  (The 
FFT-routine  automatically  adjusts  its  input  for  maximum  dynamic  range,  so 
that  the  effect  of  the  right  shift  is  essentially  to  convert  the  least  sig- 
nificant digit  to  zero.)  Ihis  was  repeated  with  a right  shift  of  2 digits, 

3 digits,  etc. 

The  test  was  applied  to  three  different  signal  sequences  derived  from 
the  USC-9  modem  with  multipath  of  a « 1/2,  t - 1.6  ms,  at  a nominal  sampling 
rate  of  7000  Hz.  Comparison  of  the  cepstra  was  restricted  to  the  region 
1 < k ^ 40  which  is  approximately  the  region  of  interest  for  multipath  detec- 
tion. Over  this  region,  the  rms  error  of  the  cepstrum  with  truncated  input 
(compared  to  the  cepstrum  without  input  truncation)  was  determined  for  each 
of  the  conditions  tried.  The  result}  are  shown  in  Figs.  B3-B6.  Figures  B3- 
B5  display  the  discrepancies  between  the  cepstrum  based  on  the  quantized 
representation  and  the  cepstrum  computed  from  the  12-bit  representation  of 
the  signal.  For  each  of  the  three  signals  and  for  any  given  number  of  binary 
digits  dropoed,  a histogram  is  presented  showing  the  number  of  cepstrum 
points  which  differed  by  0,  by  1,  by  2,  etc.  The  units  in  these  figures  are 
based  on  the  binary  logarithm;  that  is,  they  are  the  units  m which  the  Main 
Calling  Program  generates  the  cepstrum  (Sec. A. 2. 5).  The  scales  in  these 
Figures  have  to  be  multiplied  by  in2  «.693  in  order  to  bring  them  in  agree- 
ment with  the  usual  amplitude  scale  used  for  cepstra. 

This  scale  factor  was  Incorporated  in  Figure  B6,  where  the  r.m.s.  cep- 
strum error  for  each  of  the  three  signal  samples  is  plotted  as  a function  of 
the  number  of  digits  dropped  at  the  input.  While  for  the  signal  sequences 
"A"  and  "B"  there  is  relatively  little  effect  produced  by  dropping  up  to  4 
digits,  signal  "C"  shows  significant  spectrum  degradation  even  upon  dropping 
a single  digit.  One  can  assume  that  with  the  12-bit  A/D  conversion  actually 
used,  the  quantization  error  produced  little  effect  in  the  cepstra  of  signals 
"A"  and  "B",  but  a noticeable  effect  in  the  case  of  signal  "C".  Clearly,  some 
signal  sequences  are  much  more  "error-prone"  than  others. 

It  is  possible  that  signal  C (after  overlapping)  had  a sufficiently 
larger  peak  value  than  the  others,  so  that  the  FFT  routine  applied  a different 
scale  factor,  thus  essentially  introducing  additional  truncation  error.  But 
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Figure  B3.  Ceps t rum  errors  resulting  from  successively  coarser 
quantization.  Signal  sample  A. 


Figure  B4.  Cepstrum  errors  resulting  from  successively  coarser 
quantization.  Signal  sample  B. 
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1 digit 
dropped 


2 digits 
dropped 


3 digits 
dropped 


4 digits  dropped 


Figure  B5.  Cepstrum  errors  resulting  from  successively  coarser  quantization. 
Signal  Sample  C. 


it  is  more  likely  that  the  different  results  obtained  with  this  signal  are 
due  to  the  same  effect  mentioned  at  the  end  of  Section  B.l.  Thus,  ".he  DFT 
of  the  signal  C may  have  more  points  close  to  zero,  thus  aggravating  the  «• 
effect  of  slight  changes  in  input  signal. 

In  Figure  B.6  are  also  shown,  for  comparison,  estimates  of  r.m.s. 
cepstrum  error  as  obtained  from  equation  (F-2) . 

The  errors  which  are  plotted  in  Figures  B3-B5  were  in  each  case  quite 
randomly  distributed  over  the  region  of  the  cepstrum  which  wes  examined. 
The  cepstrum  locations  containing  the  multipath  spike  were  not  affected  in 
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B.4  Sampling  Jitter 

The  discussion  of  sources  of  error  contained  in  the  preceding  sections 
is  useful  in  understanding  the  overall  system  operation  and  can  serve  to 
guide  further  program  design.,  However,  in  the  experimental  setup,  the  error 
due  to  the  sources  which  have  been  discussed  so  far  was  overshadowed  by  the 
error  due  to  sampling  Jitter,  - a fact  which  was  recognized  too  late  to  be 
corrected. 

The  excessive  sampling  jitter  was  a result  of  the  manner  in  which 
synchronization  for  the  sample-and-hold  module  was  derived.  Whereas  the 
sampling  rate  was  set  by  means  of  a squarewave  generator,  the  "hold"  was 
initiated  when  a loop  in  the  A/D  subroutine  first  encountered  a "positive" 
condition  at  the  squarewave  generator  output.  Therefore,  the  initiation 
of  "hold"  had  a variable  delay  with  respect  to  the  positive-going  edge  of 
the  square  wave.  The  maximum  value  of  this  variation  in  delay  is  the  cycle 
time  of  the  software  loop  in  question,  which  is  the  following: 

Loc:  Contents  Operation:  Exec.  Time 

0307  6141  If  input  timing  * 1,  skip  next  instruction.  2.6  psec. 

0310  5307  Jump  to  0307  3 psec 

Initiation  of  "hold"  is  effected  by  the  next  instruction.  The  maximum 
variation  in  "hold"  timing  is  therefore  5.6  psec. 

The  r.m.s.  error  due  to  sampling  jitter  can  be  determined  in  the  following 

way.  Let  the  sampling  jitter  be  characterized  by  a uniformly  distributed 

random  variable,  J,  and  let  4 denote  the  peak-to-peak  value  of  the  jitter. 

If  *(t)  is  the  signal  being  sampled,  then  the  sampling  error  at  time  tQ  is 

Jx'(tQ).  (This  expression  assumes  that  x*(t)  does  not  vary  substantially 

over  a A-interval,  which  is  true  in  the  case  being  considered.)  The  r.m.s. 

Jitter  at  t is  therefore  (A/  \/T2)x’(t  ).  The  r.m.s.  value  of  the  deriva- 
o o 

tlve  of  the  input  signal  can  be  determined  approximately  from  the  signal 
sample  sequence  by  using  (B-9).  The  approximation  is  based  on  the  simplifying 
assumptions  of  Sec.  B.5.  The  approximate  m.s.  error  due  to  sampling  jitter, 
at  the  sampler  output,  is  therefore 
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(B-l) 


N - ! 

.2  * 2irnf  2 , «>  « 2 

lit2  I |dx(n)|2  + (Trf  )2|dx(|)|2). 


B. 5 Computing  the  m.s.  value  of  the  derivative  of  a signal 
from  the  sampled  signal 

Consider  a real  signal  x(t)  which  is  sampled  to  give  the  sample  sequence 

{xk}  " {x(F)}  (k  " °»  N_1>*  Its  DFT  18 


N-l  * 27rnk 


(dx(n)}  - l xke  N >• 


(B-2) 


The  "energy"  of  the  sequence  (dx(n)}  is  the  m.s.  value  of  the  signal  sequence 

N-l  N-l  N-l  2imk  N-l  +j  2irnl  N-l 

l ldx(n>|2,,i2  I ( l xke  * X l N )-N  2 Xk* 

n-0  N n-0  k-0  1-0  k-0 

CB-3) 


Equivalently, 


m.s.  value  of  DFT  - — x m.s.  value  of  signal  sequence. 


(B-4) 


The  inverse  of  the  relationship  (B-2)  is 
N-l  j 2irnk 

{x.}  m l l d (n)  e N }, 


which  can  also  be  written 


(Xfc)  " dx<°)  + I dx(n)e 
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If  the  signal  x(t)  is  assumed  periodic  with  period  T,  and  strictly  band- 
limited  to  frequencies  not  exceeding  then  x(t)  can  be  expressed  exactly 
by  replacing  kT/N  in  (B-5)  by  t: 


(B-6) 
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The  derivative  of  x(t)  is  therefore  given  by 
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The  sample  sequence  of  the  derivative,  {x'},  then  has  DFT  {d  , (n)}  given  by 


{dx,(n)>  - {j  ^dx(n)} 


, - N N. 

in  ■ x - "j,  • • • » 2'  • 


(B-8) 


Using  (B-3)  and  the  symmetry  in  the  DFT  of  a real  sequence,  the  m.s.  value 
of  the  derivative  is  found  to  be 


N-i  / f - 1 2 

l |dx,<n)|2  -|2  l (ill)  |dx(„)|2 

n-0  \ n»l  ^ 
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(B-9) 


B.6  Estimate  of  the  Mean  Signal  Censtrum 

Estimation  of  multipath  amplitude  requires  subtraction  of  the  mean  signal 
cepstrum  from  the  actually  observed  cepstrum.  Since  the  mean  signal  cepstrum 
is  only  available  ii  the  form  of  an  estimate,  it  introduces  an  error  in  the 
multipath  estimate.  This  is  a systematic  error,  for  a given  type  of  signal 
and  sampling  rate  (l.e.,  for  a set  of  multipath  estimates  all  of  which  are 
based  on  the  same  mean  cepstrum  estimate.) 

A representative  computation  of  the  magnitude  of  this  error  was  performed 
using  6 different  mean-ceps trum  estimates  which  were  obtained  with  the  USC-9 
modem  at  7000  samples/sec.  (see  Fig.  Cl).  By  averaging  over  cepstrum  points 
number  3 to  number  22  (corresponding  to  a range  of  t from  .43  ms  to  3.14  ms) 
the  average  standard  deviation  of  the  mean-cepstrum  estimates  was  tound  to 
be  1.93.  j 

Another  error  arises  due  to  the  fact  that  the  true  mean  cepstrum  is  1 
actually  modified  somewhat  if  multipath  is  present  (aside  from  the  super-  ) 
position  of  the  multipath  spikes).  This  error  increases  with  both  a and  t. 
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and  Is  due  to  the  end  effect  illustrated  in  Fig.  3.  In  order  to  estimate  the 
size  of  this  error,  it  is  necessary  to  compare  mean  cepstra  that  have  been 
obtained  with  multipath  present  (for  various  values  of  a and  t) , ard  the 
mean  cepstrum  without  multipath.  This  procedure  also  would  yield  measurements 
of  the  attenuation  experienced  by  the  multipath  cepstrum,  as  discussed  in  Sec. 
2.4.  No  systematic  sequence  of  such  tests  was  performed,  however. 

B.7  Detection  Thresholds 

The  setting  of  the  detection  thresholds  hinges  on  the  probabilities  of 
several  different  events: 

pffl,  the  probaolllcy  of  missing  the  multipath  spike  (given  that  multi- 
path  is  actually  present) 

p^,  the  probability  of  a false  alarm  (detecting  a multipath  spike  given 
that  no  multipath  is  present) 

p , the  probability  of  wrong  detection  (detecting  a multipath  spike 
w 

at  a wrong  position  in  the  cepstrum,  given  that  multipath  is 
actually  present) 

Suppose  a single  multipath  component  is  present  with  delay  t ■ rtQ  and 

relative  amplitude  a.  If  the  distribution  of  the  signal  cepstrum  at  point  r 

is  symmetrical  and  the  mean  has  been  subtracted,  then  P ■ 1/2  if  No/2  ■ 0 . 

m r 

(This  neglects  the  possibility  of  exceeding  the  reverse  polarity  threshold) . 
Thus,  as  far  as  the  probability  of  a "miss'*  is  concerned,  the  setting  of  the 
threshold  is  guided  by  how  small  a multipath  amplitude  is  to  be  detected. 

The  larger  the  true  a when  a "miss"  occurs,  the  more  serious  the  error  at  the 
output  of  a multipath  correction  network  which  has  been  adjusted  in  accordance 
with  the  detection  result. 

Taking  into  account  the  possibility  of  false  alarms  is  not  quite  so 
straightforward.  If  the  thresholds  are  set  high,  then  P^a  is  small,  but  the 
false  alarms  which  do  occur  will  result  in  significant  distortion  when  multi- 
path  cancellation  is  applied.  On  the  other  hand,  many  of  the  false  alarms 
which  occur  under  a low  threshold  setting  will  produce  small  o-estimates,  so 
that  less  serious  distortion  is  introduced  by  a multipath  correction  network 
which  has  been  adjusted  in  accordance  with  the  detection  result. 
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A wrong  detection,  o : -letectlon  error,  occurs  when  some  point  in  the 
cepstrum  exhibits  such  a large  excursion  (due  to  the  random  nature  o£  the 
signal)  that  it  is  interpreted  as  the  multipath  spike,  while  a multipath 
spike  is  actually  present  at  a different  point  in  the  cepstrum,  but  with 
smaller  magnitude.  A wrong  detection  results  in  a double  degradation  in 
place  of  multipath  cancellation  - due  to  the  uncorrected  actual  multipath 
and  due  to  the  compensation  of  a multipath  component  actually  not  present. 
Here  again,  increasing  the  threshold  setting  lowers  P . Also,'  it  allows  only 

V 

the  more  disruptive  of  such  errors  to  occur,  which  are  those  associated  with 
large  values  of  a. 

An  approximate  expression  for  is  obtained  below  by  neglecting  the 
possibility  of  one  or  more  of  the  minor  spikes  of  the  multipath  cepstrum 
falling  within  the  detection  window.  However,  the  possibility  is  included 
that  the  cepstrum  value  at  the  location  of  the  ma^or  multipath  spike  does 
not  reach  threshold. 

Suppose  that  multipath  is  present  with  a delay  equal  to  an  integer 
number  of  sample  spaclngs,  t **  r(T/N),  so  that  the  main  multipath  spike  is 
located  at  cepstrum  point  r.  Let  the  multipath  amplitude  be  a,  so  that  the 
multipath  spike  has  amplitude  c^Cr)  ■ Na/2  - 6^.  (The  parameter  $r  expresses 
the  depression  of  the  multipath  spike  at  point  r due  to  the  fact  that  multi- 
path  introduces  a linear  shift,  not  a cyclic  shift.)  The  signal  cepstrum  at 
point  r is  c^Cr).  Thus,  a detection  error  arises  if  c^(r)  ■ cm(r)  + cx(r)  < 
cx(k)  for  some  k within  the  detection  window  such  that  cx(k)  > 0^. 

The  following  simplifying  assumptions  will  now  be  made: 

(a)  c^Ck) , for  different  k within  the  detection  window,  are  independent, 

identically  distributed  random  variables;  a representative  one  of  these 
will  be  denoted  C. 


(b)  The  thresholds  for  each  k are  all  set  to  the  same  values,  +0,  with 
respect  to  the  true  mean  of  the  signal  cepstrum  at  point  k. 

Let  D be  the  number  of  cepstrum  points  lying  within  the  detection  window. 
The  conditional  probability  of  error,  given  Cy(r),  (where  c^.(r)  > ©r)  is  then 
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PCEj<^(r))  - 1 - {FcIcyCr)J  - Fc[-cyCr)}} 


CB-10) 


D-l 


The  probability  of  a detection  error,  given  c^Cr),  is  then 
00 

Pw  - f P(e|c+  cm(r))fc(c)dc  + po(l  - [Fcce)  - FcC-6)]D'1}.CB-11) 

6 -C  (r) 
r in 

In  this  expression,  pQ  is  the  probability  of  the  event  cyCr)  < 0f.  A 
slightly  simpler  expression  than  (B-ll)  is  obtained  for  the  combined  proba- 
bility of  a miss  or  a detection  error: 


P + P 
m w 


F 


er“c_Cr) 

r m 


P(E|c+  cm(r))fG(c)dc  + pQ. 


(B-12) 


If  (B-12)  is  computed  for  the  case  where  the  C's  are  gaussian  random 
variables  with  standard  deviation  o,  and  with  6 * 0,  the  following  repre- 
sentative values  are  obtained: 

0*2 a 9 * 2.5 a 


a * 3o 
D - 31 


.128 


.047 


a - 3a 
D - 16 

a * 4a 
D - 31 


.074 

.04b 


.025 

.024 


a * 4a 
D - 16 

a ■ 5a 
D - 31 

a * 5a 
D - 16 


.026 

.0072 

.0039 


.0124 

.0050 

.0026 


This  set  of  computed  values  illustrates  the  following  points: 

a)  There  is  a significant  penalty  associated  with  allowing  too  large  a 
detection  window.  The  detection  window  should  not  be  wider  than 
required  to  accommodate  the  expected  range  of  x. 

b)  The  value  of  a is  of  great  importance.  For  instance,  going  from 
{0  * 2a, a * 4a}to  {0  * 2.5a,  a * 5a)  corresponds  to  reducing  the 
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variability  of  the  signal  cepstrum  by  33Z,  and  results  in  a reduction  in 

P + P by  a factor  of  10,  approximately, 
in  w 


The  false  alarm  probability,  P^,  is  given  by  the  following  approximate 
expression,  based  on  the  above  assumptions: 


P 


fa 


1 - 


tD 


fcCc)dc 


CB-13) 
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APPENDIX  C.  RESULTS  OF  TEST  RUNS 


Cepstrum  statistics  were  collected,  and  multipath  detection  runs  per- 
formed, using  signals  from  three  different  types  of  modems:  a USC-9  modem, 
a Codex  9600  modem,  and  a Hughes  HC-276  modem.  The  test  setup  for  all 
tests  was  as  described  in  Section  A3.  Analysis  of  the  test  data  is  con- 
tained in  Appendix  D. 

C.l  The  USC-9  Modem 

The  USC-9  was  operated  at  a data  rate  (serial  binary  input)  of  2400 
bits/sec.  Its  output  consists  of  four  quadraphase  PSH  channels  centered 
at  different  frequencies  within  the  audio  band. 

C.1.1  Cepstrum  Statistics 

Cepstrum  statistics  were  collected  at  several  different  sampling  fre- 
quencies (fg  « 5800,  6000,  and  7000  Hz).  The  mean  cepstra  obtained  in  this 
manner  were  used  in  order  to  perform  multipath  detection  runs  at  each  of 
these  sampling  frequencies.  With  f ■ 7000  Hz,  six  independent  computa- 
tions of  cepstrum  statistics  were  made,  in  order  to  examine  the  variability 
of  the  estimates  which  are  obtained.  Fig.  Cl  shows  the  ranges  of  cepstrum 
maxima,  minima,  and  mean  values  for  the  six  different  computations.  The 
solid  line  connects  that  particular  set  of  cepstrum  mean  values  which  was 
used  in  the  multipath  detection  runs  at  f ■ 7000  Hz. 

C.l. 2 Multipath  Detection 

Test  runs  were  made  at  sampling  frequencies  of  f » 5800,  6000,  and 
7000  Hz.  Multipath  detection  runs  at  each  of  these  values  of  fg  utilized 
the  mean  cepstrum  computed  for  that  sampling  frequency.  The  same  threshold 
settings  were  used  for  all  detection  runs,  with  rather  low  threshold  values 
and  with  the  detection  window  extending  over  cepstrum  samples  No.  1 through 
39.  The  multipath  parameters  in  all  these  runs  were  a ■ 0.5,  x - 1.6  ms. 
The  a-estimates  which  were  obtained  are  plotted  in  Fig.  C2.  No  correction 
for  the  multipath  spike  depression  discussed  in  Sec.  2.4  was  applied  to 
this  data.  Those  runs  which  resulted  in  detection  based  on  an  incidental 
cepstrum  peak  (at  an  incorrect  value  of  x)  are  indicated  in  the  column 
marked  "detection  error". 
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Range  of  cepatrum  maxima,  minima  and  mean 
cepatra,  with  fg  * 7000  Hz. 
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Figure  C2.  Multipath  Detection  at  different  Sampling  Rates  - 

USC-9  Modem.  Estimates  of  ol,  (true  a * 0.5). 
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Figure  C3.  Multipath  Detection  at  different  sampling  rates  - 
USC-9  Modem.  Estimates  of  T,  with  etf  ■ 0.5. 
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Estimates  of  t obtained  by  app1ying  (A-5)  to  the  data,  are  plotted  in 
Fig.  C2.  These  computations  weie  done  in  APL,  since  the  PDP-8  programs  do 
not  give  t.  The  results  shown  in  Fig.  C3  can  only  be  interpreted  on  a rela- 
tive basis  for  each  sampling  rate  because  in  the  tests,  a precise  adjustment 
(or  measurement)  of  the  sampling  rate  with  respect  to  the  true  multipath  de- 
lay was  not  attempted.  Thus,  with  6000  samples  per  second  (Fig.  C3-b)  the 
multipath  delay  was  sufficiently  close  to  an  exact  multiple  of  the  sample 
spacing  that  in  none  of  the  cas ,'s  was  the  threshold  reached  at  a point  ad- 
jacent to  the  maximum  cepstrum  peak.  The  same  value  of  t was  therefore  ob- 

* A 

tained  in  all  cases,  namely  t ■ i/fg  w^ere  i is  the  position  of  the  largest 
cepstrum  peak.  (See  also  Sec.  D.l  for  further  discussion  of  the  effect  of 
the  threshold  setting  on  T.) 

C. 2 The  Codex  9600  Modem 

This  modem  uses  combined  phase  and  amplitude  modulation  of  a single  sub- 
carrier. The  modem  handles  a maximum  data  rate  of  9600  bits  per  second,  but 
it  was  operated  at  4800  bps  in  the  tests  which  were  performed.  The  variabil- 
ity of  the  signal  cepstrum  was  found  to  be  similar  as  for  the  USC-9  and 
HC-276  modems. 

C.2.1  Detection  Runs 

An  extensive  set  of  detection  runs  was  carried  out  under  the  following 

conditions:  a ■ 0.5,  t ■ 2.9  ms.,  f ■ 7000  Hz.  With  this  combination  of 

» 8 

r and  fg,  the  multipath  delay  does  not  correspond  to  an  integer  number  of 
sample  spaclngs.  The  results  are  plotted  in  Figs.  C4,  C5  and  C6.  A total  of 
201  multipath  estimates  were  computed,  of  which  25  resulted  in  a detection 
error.  'There  were  no  misses.  Fig.  C4  shows  the  distribution  of  o-estimatfes 
(including  detection  errors)  without  correction  for  multipath  cepstrum  de- 
pression. Fig.  C5  shows  the  distribution  of  T-estimates,  again  excluding 
detection  errors.  The  horizontal  scale  in  this  Figure  assumes  a sampling 
frequency  of  exactly  7000  Hz.  Under  this  assumption,  cepstrum  sample  # 2 0 
corresponds  to  t * 2.857  ms,  and  cepstrum  sample  £21  corresponds  to  i ■ 

3.00  ms.  Since  the  existing  Detection  Routine  does  not  compute  T-estlmates, 
the  data  for  Fig.  C5  was  computed  on  APL,  and  the  horizontal  scale  of  the 
plot  was  quantized  to  increments  of  .05.  The  estimates  in  the  range  from 
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Figure  C4.  Multipath  Detection,  Codex  9600  Modem. 

Estimates  of  a,  with  true  a - 0.5,  t * 2.9  ms,  f - 7000  Hz. 
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2.90  to  2.95  resulted  from  interpolation  between  two  adjacent  cepstrum 
samples  (#20  and  #21),  whereas  an  estimate  of  2.86  0.00)  ms  resulted  when 
cepstrum  sample  #21  (#20)  did  not  exceed  threshold. 

The  detection  errors  are  plotted  in  Fig.  C6,  together  with  the  detection 
thresholds  which  were  used.  These  thresholds  were  "hand-picked"  on  the  basis 
of  measured  cepstrum  statistics.  A slightly  larger  negative  threshold  could 
have  eliminated  a few  of  the  detection  errors,  but  most  of  the  errors  .are 
due  to  rather  large  random  peaks  which  can  arise  in  the  cepstrum.  Errors 
which  occur  adjacent  to  the  location  of  the  true  multipath  spike  (locations 
19  and  22)  are  in  part  due  to  the  negative  overshoot  of  the  sinc-function 
which  is  associated  with  the  multipath  cepstrum. 

C.3  The  Hughes  HC-276  Modem. 

The  HC-276  generates  a single  quadraphase  PSK  signal.  It  was  operated 
at  2400  bits/sec.  Cepstrum  statistics  were  collected  and  plotted  (Fig.  C7) 
in  the  same  manner  as  for  the  USC-9  modem.  A number  of  detection  runs  per- 
formed with  this  modem  gave  similar  results  to  those  obtained  with  the  other 
two  modems. 


APPENDIX  D.  ANALYSIS  OF  TEST  DATA 


D.l  Analysis  of  Codex  9600  Multipath  Detection  Data. 

The  Codex  9600  multipath  detection  test  at  f ■ 7000  Hz,  whose  results 
are  plotted  in  Figs.  C4,  C5,  and  C6,  is  the  only  test  in  which  a large  number 
of  measurements  were  performed.  Nevertheless,  the  analysis  given  below  is 
distorted  by  the  fact  that  detection  errors  occurred,  and  that  for  each  de- 
tection resulting  in  an  error,  it  is  not  known  what  the  a-estimate  would  have 
been  in  the  absence  of  an  error.  In  other  words,  the  analysis  is  based  only 
on  those  cr'-estimates  which  were  associated  with  r-estimates  in  the  range  from 
2.86  to  3.00  ms  (see  Fig.  C5) . 

Out  of  the  201  individual  cepstrum  computations,  25  led  to  a detection 
error  (with  the  particular  threshold  settings  used  - see  Fig.  C6) . The  re- 
maining 176  correct  detections  are  characterized  by: 

in  units  of  in  units 

cepstrum  amplitude  of  a 

average  value  58.7  .458 

standard  deviation  13.0  .102 

The  fact  that  the  average  value  is  less  than  the  true  a (a  » .5,  corre- 
sponding to  a cepstrum  amplitude  of  64.0)  is  of  no  concern.  The  discrepancy 
represents  the  multipath  cepstrum  depression  discussed  in  Section  2.4  (6^.  in 

Section  B.7).  Actually,  there  is  another  small  contribution  to  this.. effect 
due  to  the  fact  that  the  interpolation  is  not  operative  unless  a cepstrum 
sample  adjacent  to  the  maximum  exceeds  threshold.  In  Fig.  C5,  a detection 
resulting  in  x ■ 2.86  or  3.00  ms  (i.e.,  no  interpolation)  in  general  yields  a 
slightly  lower  value  of  a than  would  be  obtained  if  interpolation  had  been 
used  without  regard  to  threshold. 

If  the  cepstrum  values  at  k ■ 20  and  k * 21  had  been  recorded  for  all 
the  cepstra  computed  in  this  test  (rather  than  only  the  cepstrum  values  at 
the  position  of  the  largest  cepstrum  spike),  it  would  have  been  possible  tc 
include  the  data  from  those  cepstra  which  resulted  in  detection  errors. 

This  presumably  would  have  resulted  in  a somewhat  lower  average,  and  also  a 
slightly  larger  standard  deviation,  since  instances  of  low  cepstrum  values 
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at  locations  k ■ 20  and  k - 21  are  more  likely  to  result  in  a detection  error. 

Based  on  the  values  in  the  above  table , a comparison  can  be  made  with 
the  probabilities  of  wrong  detection  computed  in  Sec.  B.7  for  the  idealized 
Gaussian  model.  From  the  above  data  and  Fig.  C6, 

a ■ 4.5  o 

6 ■ 2 o (D-l) 

D - 31 

On  the  basis  of  the  assumptions  in  Sec.  B.7,  this  should  give  a probability 
of  wrong  detection  or  miss  of  about  .02.  Actually,  detection  errors  occurred 
in  25  out  of  201  trials,  which  is  a relative  frequency  of  .124.  Tha  discrep- 
ancy is  mainly  due  to  the  following: 

a)  The  cepstrum  variability  is  not  adequately  characterized  by  a 
Gaussian  distribution;  it  does  not  fall  off  as  rapidly.  Since  the 
detection  algorithm  searches  for  a maximum,  this  is  probably  the 
most  significant  item. 

b)  The  threshold  settings  were  not  carefully  adjusted  at  every  point 
in  the  cepstrum,  and  presumably  at  some  locations  they  were  less 
than  the  2a  in  (D-l) . 

c)  Errors  in  the  mean  cepstrum  estimates  (see  Sec.  B.6)  are  not 
accounted  for. 

D.2  Contributions  of  Different  Sources  of  Error. 

The  test  data  indicates  a rather  large  variability  of  the  signal  cepstrum, 
for  all  three  types  of  modems  used  in  the  tests,  compared  to  the  synthetic  16- 
channel  PSK  signal  which  was  the  basis  for  the  earlier  study  reported  in  [1], 
Since  the  PDP-8  routine  for  computing  cepstrum  statistics  did  not  compile 
estimates  of  the  standard  deviation  at  each  cepstrum  point,  an  accurate  com- 
parison is  not  possible.  A rough  comparison  is  possible  by  noting: 

Est.  Standard  deviation  of  256-point  16-channel  PSK  signal  (Figure  31 
of  [6]): 

Approx.  5,  over  a large  part  of  the  detection  window. 

Standard  deviation  of  multipath  amplitude  estimate,  CODEX  9600,  excluding 
wrong  detections  (Fig.  C4,  and  Sec.  D.l): 

Approximately  13. 
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The  question  is  now,  whether  the  variability  encountered  in. the  test  with  the 
Codex  modem  is  due  to  the  nature  of  the  signal  only,  or  whether  it  is  signif- 
icantly influenced  by  the  various  types  of  errors  discussed  in  Appendix  B. 

In  this  section,  the  effect  of  the  various  sources  of  error  on  the  variability 
of  the  cepstrum  will  be  estimated. 

a)  Truncation  error,  on  the  basis  of  a single  test  (Sec.  B.l),will  be 
assumed  to  introduce  an  additive  error  in  the  cepstrum  with  standard 
deviation  of  .267.  (The  units  are  such  that  a quantization  interval 
has  size  1.0.) 

b)  Signal  quantization  superimposes  an  additive  noise  at  the  input. 

This  effect  will  be  combined  with  other  noises  appearing  at  the 
input  (i.e.,  prior  to  the  nonlinearity.)  The  quantization  noise  has 
standard  deviation  1/  V 12  * 0.29. 

c)  Data  preparation  superimposes  a noise  on  the  input  (see  B.2)  with 
standard  deviation  .35. 

d)  The  error  due  to  sampling  jitter  can  be  estimated  from  (B-l) • A 
Codex  9600  signal  sample  sequence  was  entered  at  an  APL  terminal  and 
the  expression  (B-l)  computed.  The  signal  had  an  r.m.s.  value  of 
484.  Application  of  (B-l)  resulted  in  an  error  standard  deviation 
of  9.16  at  the  input  which  clearly  overshadows  the  errors  from  other 
sources.  Taking  into  account  the  effect  of  data  overlapping,  but 
neglecting  the  effect  of  the  window,  the  signal  r.m.s.  value  becomes 
484/  VT,  and  the  errors  in  b) , c) , and  d)  combine  to  give 

V(.29//5)2  + .352  + (9.16'V2)2  - 6.49  (D-2) 

Using  (F-2)  as  a rough  means  for  estimating  cepstrum  error  gives  e ■ 

c 

2.68.  This  can  be  combined  with  the  truncation  error  of  item  (a)  above, 
ausuming  that  these  errors  superimpose  in  the  cepstrum,  to  give 

V 2.682  + .2672  - 2.7  (D-3) 

Using  this  value  to  correct  the  standard  deviation  of  the  ot-estimates  for  the 
effects  of  errors  gives 

Vl3.02  - 2.72  - 12.7  (D-4) 
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This  differs  very  little  from  the  uncorrected  value  of  13.0.  However,  it 
must  be  remembered  that  (F-2)  is  not  directly  applicable  to  signals  from  the 
Codex  9600  modem;  it  applies  to  the  16-^hannel  PSK  signal  which  Tugnait  [2] 
used  in  his  analysis.  However,  use  of  this  formula  here,  without  rederiva- 
tion based  on  frequency  domain  statistics  for  the  type  of  signal  in  use,  can 
at  least  be  expected  to  yield  a "ball  park"  figure.  Also  the  comparison  in 
Fig.  B6  of  cepstrum  errors  computed  via  (F-2)  and  values  obtained  directly 
from  different  signal  samples  indicates  that  the  error  estimate  obtained  with 
(F-2)  is  low.  Fig.  B6  is  based,  of  course,  on  data  obtained  with  the  USC-9 
modem,  but  the  cepstrum  characteristics  observed  with  all  three  different 
modem  types  seem  roughly  similar.  Thus,  (F-2)  appears  to  best  approximate 
the  behavior  of  the  less  error-prone  signal  sequences.  It  can  be  concluded 
from  (D-3)  that  sampling  jitter  very  likely  contributed  noticeably  to  the 
variability  in  estimated  multipath  amplitude,  whereas  the  effect  of  other 
sources  of  error  is  fairly  insignificant. 
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APPENDIX  E.  TWO  MULTIPATH.  COMPONENTS 

E.  1 The  Cepstrum  Due  to  Two  Discrete  Multipath  Components 

The  computations  and  tests  performed  on  this  project  have  mostly, 
for  simplicity,  assumed  that  only  a single  multipath  component  is  super- 
imposed on  the  desired  signal.  At  first  examination  it  appears  that  the 
case  of  two  (or  more)  multipath  components  is  a straightforward  extension 
of  the  analysis  for  a single  multipath  component,  with  an  appropriate  modi- 
fication in  the  expansion  of  the  logarithm  (.[1],  equ.  2-12). 

Thus,  consider  a received  signal  of  the  form 

y(t)  - x(t)  + ax(t  - t1)  + gx(t  - r2).  CE-1' 

where  x(t)  is  the  transmitted,  or  desired,  component.  A sample  sequence 
{y^}  then  has  a DFT  which  is  approximately  (that  is,  neglecting  error  due 
to  end  effects  or,  equivalently,  assuming  that  the  signal  is  periodic  with 
period  T) 

(d  (n)}  - (d  (n)  (1  + ae"j27rnr/N  + Be“j2lrns/N)},  (E-2) 

y x 

where  it  is  assumed  for  convenience  that  and  t2  are  integral  multiples 
of  the  sampling  interval:  ■ r(T/N),  and  x2  “ s(T/N).  The  logarithmic 

spectrum  is  then 

Uy(n»  - <tx<n»  + Ujl  ♦ „,-.)*"«/»  + Be"J2,ms/N|}.  (E-3) 

The  pattern  of  multipath  spikes  in  the  cepstrum  arises  from  the  last  term  in 
(E-3) . Its  form  becomes  clear  upon  expanding  the  logarithm: 

£ (1  + ae’j2lTnr/N  + (Se'j2WN)  - (ae_J2nnr/N  + 0e~J  2WN) 
n 

- | (a2e"jW/N  + 2aee'J2irn(r+s)/N  + 0Vj4ffns/N) 

+ (E-4) 

The  following  cepstrum  components  are  therefore  obtained: 
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(E-5) 


/ 


/ 

/ 


However,  the  expansion  (E-4)  only  converges  for 

1*1  + lei  < i, 

and  it  was  first  shown  in  [7]  that  the  amplitudes  of  the  cepstrum  spikes 
have  a completely  different  behavior,  as  functions  of  a and  g,  once  the 
condition  (E-5)  is  exceeded.  In  that  case,  the  amplitude  of  cepstrum 
spikes  can  be  obtained  numerically  from  cepstrum  computations  as  shown  for 
instance  in  Fig.  El,  which  is  a reproduction  of  Fig.  5 in  [7].  Here  the 
coefficients  of  significant  cepstrum  components  have  been  plotted  as  a func- 
tion of  a for  the  case  a ■ g,  with  r ■ 10,  s ■ 17.  For  a ■ g < 0.5,  the 
linear  increase  in  fundamental  components  and  quadratic  increase  in  second- 
order  sum  components  can  be  observed,  in  agreement  with  the  above  chart.  A 
drastic  change  occurs  in  the  pattern  near  a ■ g ■ 0.5.  The  fundamental  and 
second-order  sum  components  become  depressed  as  a and  g increase  further, 
and  instead,  difference  components  arise  at  s - r 1 7,  2r  - s * 3,  2s  - r “ 

24,  etc.  The  same  phenomenon  is  also  illustrated  in  Fig.  E2,  which  is  a 
plot  similar  to  the  one  in  Fig.  El,  but  for  the  case  where  0 ranges  from  0 
to  1 while  a is  fixed  at  a * 0.8.  As  0 increases  above  0.2,  (E-5)  is  violated, 
the  cepstrum  spikes  at  r,  2r,  and  3r  (cepstrum  points  10,  20  and  30)  get  de- 
pressed and  the  difference  component  at  s - r (*7)  appears.  Figures  El  and 
E2  show  only  the  major  cepstrum  components  which  appear  in  the  portion  of 
the  cepstrum  that  is  illustrated.  Additional  weak  components  arise  at 
higher  order  differences  when  (E-5)  is  violated. 

E. 2 Estimation  of  Multipath  Parameters 

A preliminary  version  of  an  algorithm  for  identifying  two  multipath  com- 
ponents was  developed  and  is  presented  below.  A PDP-8  program  implementing 
the  algorithm  has  also  , written,  but  not  tested.  A discussion  of  the 
algorith;  follows  its  description. 

E.2.1  Description  of  Algorit'^i 

Part  1:  Multipath  Identification 

a)  Identify  largest  cepstrum  spike  which  exceeds  the  threshold,  and  interpo- 
late. Denote  its  interpolated  peak  value  A^  and  its  interpolated  cepstrum 
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location  A^.  If  there  Is  no  point  at  which  threshold  is  exceeded,  the 
decision  is  "No  Multipath"  and  the  routine  terminates.  Otherwise,  con- 
tinue to  step  b. 

b)  Identify  the  next  largest  cepstrum  spike  which  exceeds  threshold,  and 

interpolate.  Denote  the  Interpolated  peak  value  A2  and  its  Interpolated 

cepstrum  location  If  no  such  cepstrum  spike  is  present,  the  decision 

is  "single  multipath",  anJ  the  multipath  parameter  estimates  are 

a - 2A./N,  T -A,/f  . Otherwise,  continue  to  step  c. 

" 1 • Is 

c)  Identify  the  third  largest  ctpatrum  spike  which  exceeds  threshold,  and 
interpolate.  Denote  the  interpolated  peak  value  A^*  and  its  interpolated 
cepstrum  location  Ly  If  no  such  cepstrum  spike  is  present,  go  to  step 
g.  Otherwise,  continue  to  step  d. 

d)  Order  A,,  A„,  A,  in  ascending  order  and  denote  them  A , A.  , A . (Thus, 

X 2 j a b c 

A < A.  < A .)  Let  6 denote  the  time  between  samples.  If 
a — b — c 

Afl  + A.  ■ A + d,  go  to  step  f.  Otherwise,  go  to  step  e. 

e)  If  A ^ 2A,  + 6,  and  also  A f*  2A  +6,  discard  A..  and  A,  and  go  to  step 

c b — c a J J 

g.  Otherwise  set  ■ Aa»  *2  “ compute  the  cepstrum  value  at  Tj  ” T j 
and  go  to  step  h. 

f)  If  A^'Aj'A^  > 0,  set  ■ A^,  t2  “ Ac«  take  A^  as  the.  cepstrum  value  at 
Tj  “ T^,  and  go  to  step  h.  Otherwise,  discard  A£  and  set  ■ Ag,  Tj  ■ 
A^»  compute  the  cepstrum  value  at  T2  ~ t^,  and  go  to  step  h. 

g)  Denote  Ag  - maxtA^.Aj)*  A^  • min(.}^,A2)*  If  A^  - 2Afl  + d/2,  and  A^  « 

9 

-A  /N,  then"single  multipath"  is  decided,  with  t » A , a * 2A  /N.  Other 
wise,  continue  to  step  h. 

Part. 2:  Double  Multipath  Parameter  Estimation 

h)  To  determine  whether  or  not  (E-5)  is  satisfied,  compute  the  value  of  the 
cepstrum  at  12“  T ^ . Denote  the  value  A^.  If  A^  exceeds  a preset 
threshold  level,  the  decision  is  (E-5)  is  violated.  In  that  case,  go  to 
step  i;  otherwise  go  to  step  j. 

i)  "Region  I"  estimates  (i.e.,  |a|  + |g|  < 1): 

«1  " Aa’  -2  " Ab’  * " 2Afl/N,  § - 2Afe/N 
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J)  "Region  II"  estimates  (i.e.,  |a|  + 1 0 1 > 1): 

Simultaneously  solve  the  following  two  equations  for  a and  (3: 

44a  - 7 .'8  - A - 83.5 
a 

448  - 77a  - Afe  - 83.5 

Retaining  a/g,  scale  a and  g to  satisfy 

a + 8 - A./57  + 1. 

~ ~ A 


Figure  E3.  Number  of  Cepstrum  Spikes  which  have  Amplitude  greater 
than  25,  for  different  values  of  a and  8.  (256-point 
real  cepstrum.) 
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E.2.2  Discussion  of  Algorithm 


The  algorithm  presented  In  the  preceding  Section  consists  of  tvo  parts. 
The  first  part  serves  to  identify  the  type  of  multipath  structure  which  is 
present,  and  to  identify  the  pertinent  cepstrum  components.  The  second  part 
provides  the  actual  pat  aeter  estimates  (in  the  case  of  tvo  multipath  compo- 
nents) . 

Steps  a and  b are  straightforward.  Interpolated  peak  values  and  loca- 
tions are  obtained  in  the  manner  described  in  Sections  A. 2. 6.1  and  A. 2. 6. 2. 
In  step  a it  is  decided  whether  or  not  multipath  is  present.  In  step  b,  if 
it  is  decided  that  only  a single  multipath  component  is  present,  the  multi- 
path  parameters  are  estimated  in  the  same  way  as  in  Section  A. 2. 6. 2. 

In  step  c,  a third  cepstrum  spike  is  sought  which  might  provide  addi  - 

tional  information.  If  a third  cepstrum  spike  is  not  found,  the  possibil- 
ity must  be  considered  that  the  tvo  spikes  which  have  been  found  are  the 
first  two  spikes  of  a single  multipath  cepstrum  (Fig.  1);  this  is  done  in 
step  g. 

Steps  d through  f apply  in  the  case  where  three  multipath  spikes  are 
detected.  It  is  then  necessary  to  decide  which  three  they  are  (see  Fig.  E3) . 
This  problem  exists  even  in  the  absence  of  variability  due  to  the  signal 
and  due  to  additive  noise,  and  it  becomes  more  difficult  when  such  varia- 
bility is  taken  into  account.  In  step  d,  if  the  relationship  A^  + A^  ■ 

A + 6 is  satisfied,  it  could  mean  either  that  the  spik,  i of  location  A and 

A^  are  the"main"  multipath  spikes,  and  A£  is  the  sum  spike;  or  that  the 

main  spikes  are  at  A.  and  A ,and A is  the  loaation  of  the  first  difference 

oca 

spike.  A decision  between  these  tvo  possibilities  is  reached  in  step  f. 

In. step  e,  the  possibility  that  the  spike  at  Ac  is  a second  harmonic  compo- 
nent is  weighed.  If  it  is  decided  that  this  is  the  case,  this  component  is 
discarded  because  it  is  not  used  in  the  parameter  estimation  procedure. 

It  should  be  noted  that  certain  combinations  which  might  arise  have 
not  been  accounted  for  in  this  algorithm  (they  could  be  included  without 
excessive  effort),  notably,  the  possibility  of  the  second  main  spike  coin- 
ciding with  the  harmonic  of  the  first  multipath  component  - l.e.,  the  tvo 
multipath  delays  having  a ratio  of  2:1. 
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The  second  part  of  the  algorithm  requires  first  a decision  as  to 
whether  (E-5)  is  satisfied  ("Region  I"  of  the  a, 6-plane)  or  not  '‘'Region 
II"  of  the  ct, 6-plane).  Unfortunately,  it  is  not  possible  to  make  a very 
sharp  discrimination  between  the  two  regions,  as  is  evident  from  Figs.  El 
and  E2.  The  decision  is  based  on  the  value  of  the  difference  component 
(which  ideally  is  zero  in  Region  I).  The  Region  I estimate,  in  step  i,  is 
straightforward  since  (E-4)  applies. 

The  procedure  for  obtaining  the  Region  II  estimate  is  based  on  detailed 
plots  of  the  multipath  cepstrum  amplitudes.  From  these,  linear  approxima- 
tions to  the  empirical  data  were  constructed  which  resulted  in  the  equations 
in  step  j.  The  last  part  of  step  j is  an  attempt  to  improve  the  estimate 
near  the  boundary  of  Region  II. 
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APPENDIX  F.  EFFECTS  OF  ADDITIVE  NOISE 


The  effect  in  the  cepstrum  due  to  additive  noise  in  the  time  domain  was 
considered  by  Tugnait  [2].  He  proposed  a simplified  mathematical  model  for 
estimating  the  r.m.s.  error  in  the  cepstrum  due  to  several  types  of  additive 
noise,  namely,  impulsive  noise,  white  Gaussian  noise,  and  quantization  noise. 
He  compared  results  obtained  with  this  model  against  tests  performed  with 
synthesized  signal  and  noise  sequences  on  APL,  and  found  reasonably  good 
agreement.  The  various  types  of  noise  which  were  considered  have  in  common 
the  property  that  their  spectrum  is  uniform  (in  a statistical  sense);  and  it 
was  found  that  the  error  introduced  in  the  cepstrum  was  similar  for  the  three 
types  of  noise,  for  given  input  signal-to-noise  ratio.  Of  course,  the  nature 
of  the  interference  produced  by  additive  noise  depends  strongly  on  the  proba- 
bility distributions  governing  the  magnitudes  |d^(n) | of  the  DFT  of  the  signal 
sample  sequences  (x^)*  Tugnait  used  in  his  analysis  and  computations  the  16- 
channel  PSK  signal,  at  8175  sample/sec,  which  was  also  used  in  the  work  de- 
scribed in  [1]. 

Tugnait's  simplified  mathematical  model  is  obtained  by  partitioning  the 
discrete  frequency  domain  into  two  sets  of  points,  one  consisting  of  those 
points  where  |d^(n) | has  a large  expected  value,  and  the  other  consisting  of 
the  rer:  ining  points.  Within  each  of  these  two  regions,  the  distribution  of 
|d^(n) | i«  rhen  assumed  independent  of  n and  is  approximated  by  a mathemati- 
cally simp.-  • -ion  (triangular).  After  utilizing  another  simplifying 
assumption,  namely  hat  the  first  ^ DFT  magnitudes  are  characterized  by 
mutually  independent  rand ’r  variables,  Tugnait  obtains  expressions  from 
which  estimates  of  the  r.m.s.  error  in  the  cepstrum  domain  can  be  computed. 

A simple  empirical  expression  can  be  fitted  to  the  error) data  obtained 
by  Tug:„  . case  of  an  interfering  impulse  (Fig.  4.6a  and  4.6b  in  [2]) 
as  follows:  1 

- 4 In  (1  + 32  Pi)  (F-l) 

“ ratio  of  r.m.s.  noise  to  r.m.s.  signal  at  the  input  > 

(averaged  over  the  input  sequence)  1 

■ r.m.s.  error  per  point  in  the  cepstrum  (averaged  over  the  whole 
cepstrum) 
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where 


/ 


Although  this  expression  is  obtained  for  the  case  of  a sirgle  interfering 
impulse,  it  appears  from  other  computations  presented  in  [2]  that  white 
Gaussian  noise  produces  a similar  cepstrum  error,  for  a given  input  noise-to- 
signal  ratio.  Therefore,  (F-l)  is  roughly  applicable  for  any  noise  with 
constant  spectrum  superimposed  on  a signal  sequence  of  the  type  considered  by 
Tugnait.  If  the  noise  is,  in  fact,  a single  interfering  impulse,  can  alsc 
be  written  ^f(S/A),  where  S is  the  signal  r.m.s.  value,  A the  impulse  magni- 
tude, and  N the  number  of  points  being  processed. 


The  expression  (F-l)  is  specifically  for  128-point  signal  sequences. 
However,  it  can  be  modified  to  apply  for  arbitrary  N in  the  following  manner. 
Consider  doubling  N,  while  keeping  p^  constant.  The  signal-to-noise  ratio  in 
the  DFT  domain  (and  log  spectrum)  is  unchanged.  (If  the  noise  is  impulsive, 
this  assumes  the  number  of  observed  Impulses  doubles.)  However,  in  computing 
the  cepstrum,  integration  (actually  summation)  is  now  over  twice  as  many 
points  of  the  log  spectrum  as  before,  resulting  in  an  improvement  in  the 
ratio  of  signal  cepstrum  to  noise  cepstrum  by  a factor  of  'fl.  Since  the 
signal  cepstrum  increases  by  a factor  of  2,  ec  increases  only  by  /2.  The 
appropriate  modification  of  (F-l)  is  therefore 


c 


c 


—7=  ln(l  + 32  p.) 
2V2  1 


(F-2) 


This  equation  is  solved  below  for  several  different  values  of  p^  with  N ■ 256: 


c 

c 


0.0001 

0.001 

0.01 

0.1 

1 


.0181 

.178 

1.57 

8.1 

19.8 


APPENDIX  G.  SSB  AND  DSB  S GNALS  SUBJECT!  TO  MULTIPATH 

G.l  Ordinary  Synchronous  Detection  of  SSB  Signal 

Suppose  the  baseband  signal  x(t)  is  used  to  SSB  modulate  a carrier  at 
frequency  fQ.  Let  9^(f)  and  ^(f)  denote  the  positive  and  negative  fre- 
quency portions  of  the  F.T.  of  x(t) , respectively,  and  xCt)  the  Hilbert  trans- 
form of  x(t).  Assuming  that  the  modulation  resides  In  the  upper  sideband, 
the  modulated  carrier  is  then  of  the  form 


v(t)  ■ x(t)  cos  w t - x(t)  sin  u t, 
o o 


where 


9\f(f)  - ?t(f  - f ) + Tl( f + f ) 


(G-l) 


(G-2) 


Introduction  of  a single,  fixed,  discrete  multipath  component  of  relative 
magnitude  a and  delay  t results  in  the  distorted  SSB  signal 

w(t)  - v(t)  + av(t  - x)  (G-3) 

whose  F.T.  Is 

T(f)  - Tv(f)  + o 7v(f)e"j2ltfT 

- T*(f  - fo)[l  + oe'j2nfT]  + 7^(f  + fo)[l  + ae“J2nfT].  (G-4) 

Assume  that  the  demodulator  is  synchronized  to  the  phase  of  a reference 
carrier  included  In  the  SSB  signal.  If  the  received  carrier  phase  angle  in 
the  absence  of  multipath  is  taken  as  0°,  then  the  actual  received  carrier 

(G-5) 


phase  is  -ju0T 

6 - arg[l  + ae  ] » arctanryy 


-a  sin  u0T 


a cos  (oqt 


Neglecting  constant  multipliers,  the  demodulated  signal  y(t)  is  therefore 
given  In  the  frequency  domain  by 

Ty(t)  - i*rw(f)*[ej05(f  - fQ>  + e‘j06(f  + fQ)]>  H^f) 

- <ej0  Tw(f  - f0)  + e'J0Tw(f  + fo)>  H^tf) 

iJMf+f  )x  _ -J2ir(f-f  )x  * 

- f^(f)[l  + ae  ° ]e  j0  + T (f)[l  + ae  ° )^0 

x (G-e: 
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here  H^Cf)  is  a low-pass  filter  function  and  6 is  given  by  (G-5) . 
ubstitution  of  (G-5)  into  (G-6)  leads  to 


(O  ■ 7 (f)(  Icos  0 - jsgnfsinB]  + {ae  ^ir^T(cos(u  x+0) 
y x o 


- jsgnf  sin(a»oT  + 0))]}, 

giving 

y(t)  ■ x(t)cos  0 + x(t)sin  0 + ax(t  - t)cos(u>ot  + 0) 

+ ax(t  - T)sin(u  t+0). 

o 

This  can  be  written 

y(t)  * t*(t)  + o*x(t  - t) ] cos  0 + [x(t)  - x(t  - t)]sin  0, 


where 

a'  ■ a(cos  hi  t - sin  w x tan  0) 
o o 


a(ci  + cos  to  t) 
o 

1 + f 108  U T 
O 


(G-7) 


(G-8) 


(G-9) 

(G-10) 


From  (G-8)  or  (G-9)  it  can  be  seen  that  this  type  of  SSB  demodulation,  in  the 
presence  of  simple  multipath  at  r.f.,  results  in  both  in-phase  and  quadrature 
components  of  the  baseband  signal.  The  in-phase  and  quadrature  delay  compo- 
nents appear  with  different  relative  magnitudes.  Also,  the  phase  angle  0 
defined  in  (G-5)  is  very  sensitive  to  small  changes  in  x,  and  so  is  a'. 

These  features  make  cepstrum  processing  for  multipath  removal  not  directly 
applicable  to  the  demodulated  signal  y(t)  obtained  by  ordinary  synchronous 
detection. 


G. 2 Complex  Demodulation. 

Complex  demodulation  of  an  SSB  signal  is  achieved  by  means  of  complex 

sampling  of  the  IF  output,  or  by  using  both  an  in-phase  and  a quadrature 

carrier  to  perform  two  synchronous  demodulations.  With  w(t)  as  defined  in 

(G~3),  let  w(t)  denote  the  analytic  signal  corresponding  to  w(t).  Neglecting 

constant  multipliers,  the  in-phase  component  of  the  output  is  then 
~ -J2wf  t 

y.  ■ $£,[w(t)e  ° e ^ ] - w cos(w  t + 0)  + w sin(u  t + 0).  (G-ll) 

1 o o 

Similarly,  the  quadrature  component  of  the  output  is 
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Let 


Ia[w(t)e 


-je,  - 

! J " W 


cos(w  t + 8)  - v 


sin(to  t + 
o 


e) 


y denote  the  complex  envelope  of  w. 

O **  “j2TTf  t m _i«  « “ j (fl+0)  T ) 

y ■ y*  + jy  ■ we  ° e 1 - x(t)e  + ax(t  - x)e 

1 q 

- [x(t)  + a'x(t  - T)]e-^6  , 


where 


-ju)oT 

ae 


(G-12) 


(G-13) 

(G-14) 


Thus,  the  complex  envelope  of  the  complex  demodulator  output  retains  the 

general  form  of  a signal  disturbed  by  a single  multipath  component.  However, 
o 

y is  a complex  function,  so  that  cepstrum  processing  is  needed  for  both  y^ 
and  y^.  And  not  only  the  polarity,  but  the  phase  angle  associated  with  o' 
must  be  determined. 


G.3  Synchronous  Demodulation  of  DSB  Transmission. 

In  the  case  of  synchronous  demodulation  of  DSB 

transmission  we  have. 

analogous  to  (G-l)  through  (G-4): 

v(t)  - x(t)  coswot. 

(G-15) 

%«>  - i + j y * + «„> 

(G-16) 

ru(f)  -rvcnn  + ae~J27rfti 

(G-17) 

The  received  reference  carrier  phase  is  again  given  by  (G-5) , 3nd  we  have 

yf>  - { *j\«  - + *'JI * *„>}  HLP(f> 

”j2ir(f+f  )t  _ 


" 2 rx<£>[1  + oe 


1 -10  1 “j  2ir  (f-f  )t 

Je  J0  +f  ^(f)[l  + oe  ° }eJ0 


- T(f)[cose  + ae“J2irfT  cosCe  + a,  T>  1 
A o 

, 2 

or  Ct  Sin  (|)  x 

■ tr  (f)  cosell  + ae  ' nfT(cos  u t + T~Z 2 33 

x o 1 + Q cos  U)  x 


" ^(f)cose{l  + ae"J2irfx 


cos  o>0t  + a 
1 + a cos 


] 


CG-18) 
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Thus,  the  received  signal  is  of  the  form 


y (t)  = x(t)  + a'x(t  - t),  CG-19) 

where  a'  is  again  given  by  (5-10) . Thus,  for  a given  a,  o'  lies  in  the 
range 

— i a | <.  a'  <.  |a  |,  CG-20) 

its  exact  value  determined  by  u t.  The  coefficient  cos  0 in  CG-18)  can  also 

o 

be  written 


cos  0 « 


1 + acosrn  t 
o 


V 1 + 2acosw  t + 


(G-21) 


G.4  Complex  Demodulation  of  DSB  Transmission. 

With  w(t)  defined  as  the  inverse  transform  of  (G-17),  let  w(t)  denote 

the  analytic  signal  corresponding  to  w(t): 
ju  t jtu  (t-x) 

w(t)  ■ x(t)e  + ax(t-x)e  (G-22) 

Again  neglecting  constant  multipliers,  the  in-phase  and  quadrature  components 
of  the  output  are  given  by  (G-ll)  and  (G-12) , respectively.  The  complex 
envelope  of  w(t)  is 

-j(u  t+0)  , 

y(t)  - y±  + jyq  ■=  w(t)e  ° - [x(t)  + o’xCt-Tjle"30,  (G-23) 

where  a'  ■ ae  , as  in  (G-14). 

As  with  complex  demodulation  of  SSB,  a*  is  complex.  However,  the  trans- 
mitted signal  appears  as  the  real  signal  in  (G-23) , whereas  it  appears  in  its 
analytic  form  iu  (G-10). 


69 


/ 


REFERENCES 


1.  H.  S chwarz lander , "Use  of  the  Cepstrum  for  Processing  Multipath 

Signals"  Technical  Report  RADC-TR-72-44 , Rome  Air  Development 
Center,  Griffiss  Air  Force  Base,  N.Y.  March  1972.  AD-740586. 

2.  J.K.  Tugnait,  "Effects  of  Additive  Noise  on  Cepstral  Analysis  of 
Signals,"  Technical  Report  TR-73-6,  Dept,  of  Electrical  and  Computer 
Engineering,  Syracuse  University,  Syracuse,  N.Y.  June  1973. 

3.  W.N.  Sirgany,  "The  Design  of  Signals  for  Effective  Filtering  by 
Cepstrum  Techniques,  with  Application  to  Multipath  hannels"  Ph.D. 
Dissertation,  Dept,  of  Electrical  and  Computer  Engineering,  Syracuse 
University,  Syracuse,  N.Y.  January  1974.  (To  be  published  as  RADC 
Technical  Report) . 

4.  R.W.  Schafer,  "Echo  Removal  by  Discrete  Generalized  Linear  Filtering" 
Ph.D.  Dissertation,  M.I.T. , Cambridge,  Mass  , February  1968. 

5.  James  Rothman,  "FFTS-C  — - A Fast  Fourier  Transform  Subroutine  for 
Complex  Bata,"  DECUS  Program  Library  No.  8-144.  August  7,  1968. 

6.  C.  Schenker,  "The  Effect  of  Signal  Sequence  Size  on  the  Ability  of 
a Cepstrum  Signal  Processor  to  Detect  Multipath"  M.S.  Thesis, 

Syracuse  University,  May  1972. 

7.  H.  Schwarz lander,  "A  Cepstrum  Processor  for  Signals  with  Multipath 
Distortion,"  Signal  Processing  Conference,  Erlangen,  Germany, 

April  4-6,  1973. 

8.  W.  Sirgany,  "Signal  Design  for  Improved  Multipath  Identification  in 
the  Cepstrum"  Tech.  Memo.  TM-73-2,  Department  of  Electrical  and 
Computer  Engineering,  Syracuse  University,  Syracuse,  N.Y.  June  1973. 


70 


s»py^,w,;| 


MISSION 

of 

Rome  Air  Development  Cente. 

RA DC  is  the  principal  AFX  organization  charged  with 
planning  and  executing  the  USAF  exploratory  and  advanced 
development  programs  for  information  sciences,  intelli- 
gence, command,  control  and  communications  technology, 
products  and  services  oriented  to  the  needs  of  the  USAF. 
Primary  MX  mission  areas  are  communications,  electro- 
mag.ietlc  guidance  and  control,  surveillance  of  ground 
and  aerospace  objects,  intelligence  data  collection  and 
handling,  information  system  technology,  and  electronic 
reliability,  maintainability  and  compatibility . MX 
baa  xdssion  responsibility  as  assigned  by  AFX  for  de- 
monstration and  acquisition  of  selected  subsystems  and 
systems  in  the  intelligence,  mapping,  charting,  command, 
control  and  communications  areas. 


